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I. INTRODlKrPIOH 


Sound is both a ptQrslcal phenomenon defined as the wave oscilla- 
tion of a medium and a physiological phenomenon defined as a sensation 
received by the esur due to acoustic waves. As a mechanical vibration 
sound is characterized by frequency, anqplitude, and phase angle. How- 
ever, no simple and unique relationship exists between the physical 
meas\irement of sound and the human perception of sound. Consequently, 
many physiological methods exist to evaluate various aspects of the 
perception of soxuid. For exanqple, the human perception of piu*e tones 
is described by normal equal loudness contours. 

Furthermore, sound is edso a psychological phenomenon which can be 
characterized in terms of emnoyance and desirability. An laidesirable, 
unwanted, or excessive sound is called a noise. Besides producing 
annoyance, noise can produce various degrees of permanent hearing loss. 
Also, since people prefer to live in quiet places and to be near quiet 
objects, noise ceui have an economic importemce. 

Due to the increasing public concern with the environment and to 
the growth of commercial aviation in this country, aircraft noise is 
considered to be one of the more significeutt noise sources. As a re- 
sult ed-rcraft noise abatement is a major objective of government and 
industry. 

Airport community noise can be reduced by applying restrictive 
aircraft and airport operational procedures or by developing new tech- 
nology. Safety and economic considerations reduce the effectiveness of 
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such operational restrictions as the nulti-segment approach, power cut- 
back, runway restrictions, establishment of airport buffer zones, and 
land use control. Consequently, in recent years noise has become an 
important consideration in the de.>ign of t\u‘bofan aircraft engines. 

In the study of turbofan aircraft engine noise, the noise is con- 
sidered to be d\ie to various sources and each source is aissociated with 
specific components of the engine or specific air flow conditions. The 
noise firom turbofan engines is often attributed to the following three 
sources; the fan, the jet exhaust , and the engine core. Fan noise prop- 
agates to the ground from both the engine inlet and exit duct planes. 
Also, Jet mixing noise is produced by both the fan and core exhaust air. 
These turbofan engine noise sources are shown schematicilly in Fig- 
iire 1. The relative magnitude of fan noise 8«id Jet exhaust noise is 
generally greater than core noise for most operating conditions of cur- 
rent engines. 

Fan noise is characterized by a high pitched whine produced by the 
rotating fan blades. The high velocity Jet exhaust generates noise as 
it mixes with the relatively quiescent atmosphere. The term "core” 
noise refers to all the other noise sources including the following: 

1. Combustion noise which is the noise associated with pressure fluc- 

tuations due to uneven bui*ning 

2. Entropy noise which is the noise due to nonuniformities in the ex- 

panding gases as they pass through the turbine stages and the 
exhaust nozzle 

3. Rotating machinery noise due to the compressor and turbine 

^4. Flow noise due to air passing through ducts and over internal sur- 


faces. 
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Nost of the noise from current turbofan engines is attributed to 
the fan or Jet. Over the past few years aircraft system quieting tech* 
nology has been developed which reduces the noise fr<» these two 
soxirces. However, research on core noise and combustion noise, in par- 
ticular, has not been neglected. 

Combustion noise research is conducted because it could be iq^r- 
tamt in future aircraft design for three reasons. First, it coul<^ be 
the dcoiinant noise source in the duct-burning turbofan engine, which is 
one of the proposed commercial supersonic aircraft engine propulsion 
systems (Stewart and Weber (1978)). Second, as the noise from the fan 
and Jet sources is reduced on future aircraft, the in^rtance of such 
secondary sources associated with the engine core as combustion noise 
increases. Third, core noise including combustion noise can be a sig- 
nificant contributor to overall engine noise during takeoff or approach 
when the fan noise and the Jet noise are reduced because of forward 
velocity effects (Stone (1975)). 

In order to characterize aircraft engine combustion noise the fol- 
lowing items are of interest: 

1. The acoustic power level variation with burner gecanetry and per- 

formance parameters 

2. Tlie transmission loss across the turbine and at the engine exit 

3. The combustion noise directivity pattern of far field engine data 
1*. The combustion noise spectral shape 

5. The peak frequency of the combustion noise spectrum. 

These are also the items a combustion noise prediction scheme would 
yield as part of an aircraft engine design procedure. 
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In ortksr to obtuln thin irifomation* Many different types of test 
fecilitiee end & veriety of MMiirenent tecluiiquet ere used to ett^ 
ecathuetion noise. Iliese tests ere conducted vlth the eiqpectetiw thet 

(1) c^bMtion noise Beesureaents on non-engine test feellities ere re- 
leted to cottbustlon noise saesuresMnts on engine test feellities end 

(2) that c«Bbu8tioD noise swasuresMnts snide using different engines ere 
related to one another. Soae of these test facilities end saesureaent 
techniques vill be described next. 

The studies of Bonnell> Marshall, and Riecke (]97l)* Mathews end 
Rekos (1976), Karchaer and Reshotko (1976), ShivasNmkara (1977)* 
Reshotko, Karchaer, Penko, and McArdle (1977), and Kercltter (1978) were 
made using operating turbofan engines. 

Combustors installed in open flow ducts were used by Katin end 
EBBoerling (197^*), Strahle, Muthukrishnan, and Heale (l977e, 1977b), 
Abdelhaaid, Harrje, Plett, and Summerfield (1973), and Plett, Leshner, 
and Sunaerfield (1975). 

The studies of Essnerling (1975), Eawerling and Bekofske (1976), 
Sofrin and Ross (1975), Sofrin and Riloff (1976), and Reshotko and 
Karchmer (1977) ’-.ere made using combustors installed in closed duct sys- 
tems called component development facilities. 

Combustion noise from open flames was measured by Smith and Kelham 
(1963), Hurle, Price, Sugden, and Thomas (1968), Knott (1971), and 
Muthukrishnan, Strahle, and Handley (1976). 

External and internal pressure measurements can be made on tests 
of combustors mounted in turbofan engines and in open flow ducts. Only 
internal pressure measurements can be made in tests of cc^ustors in- 
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•tfdled in closed d\»ts» and (mly external pressure »MSursMiits can be 
BMtde of the noise of open flaswe. 

In the study of eomhfustion noise, pressure ■easursnents are sub- 
ject to one-third octave band analysis or narrow band spectrua analysie 
to obtain frequency and asq>litude infonsation. Also, transfer function, 
correlation, and coherence calculations are s«iMtlBes Bade between pairs 
of pressure ateasuresents to obtain inforoation on Mplitude and i^ase 
changes dttring propagation, on time delays between signals and on the 
signal to noise ratio. 

Experimental studies of combustion noise using operating turbofan 
engines are made with re^istic ccssbustion region boundary conditions 
and at reaJistic ccnabustion region operating conditions. However, sep- 
arating combustion, fan, compressor, turbine, and Jet noise contribu- 
tions from one another in external far-field acoustic measurements can 
be difficult. Consequently, correlation and coherence techniques using 
internal and external microphones have been developed to identify co»- 
bustion noise in the far-field acoustic measurements and to quantify its 
contribution (Karchra'^r and Reshotko (1976), Shivashankara (1977), 
Reshotko, Karchmer, Penko, and McArdle (1977), Karchmer, Reshotko, and 
Montegani (1977), Karchmer (1976)). Using these techniques, the com- 
bustor has been shown to be a source of low frequency core noise 
(Karchmer and Reshotko (1976), Karchmer, Reshotko, and Montegani 
(1977), Karchmer ( 1976 )). 

Extensive testing of many combustor desi^s on full-scale turbofan 
engines using correlation and cohere, e techniques can be expensive. An 
alternative procedure is to conduct acoustic design studies on a less 
expensive teat facility. The combustor component duct rig is one type 
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of less expensive facility. Because this facility is closed, operating 
the combustor at realistic conditions of teopervture* pressure, and flow 
is possible. However, the conbustion region boundary conditions pro* 
duced by a cosibustor ensponent duct facility differ from those of an 
engine since the turbine is not included in the easQ)onent dwt rig. 

This luy not be too significant since von Glahn (1978) was able to cor* 
relate combustor acoustic power levels inferred from internal fluctua* 
ting pressure neasurenents made on both the engine and the combustor 
cosponent duct facilities with operating conditions and chamber geome* 
tries. 

Constructing a combustor test facility which costs less than an 
engine can also be accomplished by installing a combustor in an open 
flow duct. For this facility also, the combustor boundary conditi».ns 
differ from those of an engine. In addition, while the combustor in 
this type of facility can operate at turbofan engine combxistor temper- 
ature and flow conditions, the combustor operating pressure for a given 
duct exit termination (horn, open duct, flange, nozzle) is also deter- 
mined by the flow and temperature specified. Consequently, matching 
engine operating conditions with this rig is difficult. Therefore, the 
open flow duct rig is used more often in conducting such fundamental 
research as separating hydrodynamic entropy, and direct combustion 
noise contributions to internal fluctuating pressure measurement 
(Strahle, Muthukrishnan, and 'Jeale (1977)) than in making parametvic 
studies. 

This disseration describes an investi^tion of the propagation of 
combustion noise in an open flow duct. This is a new area of re- 
search. The research is 


conducted using a combined experimental and 
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theoretical approach, l^ie objective of the theoretical investigation is 
the study of the interaction of the coabustion noise propagating in a 
plane wave stode with the duct combustion environment. In this case* the 
duct combustion environment includes the emissions produced by the cosh 
bustor* the air flow, and the duct exit impedance and inlet isgtedance. 
Ihe interaction of a propagating plane acoustic wave with eoabxistor 
emissicns has not previously been considered. This work uses a dynamic 
systems state-space approach to study this problem. 

This investigation was initiated due to difficulties encountered 
in the study of the structure of the internal fluctuating pressure 
spectra measurements taken in an open flow combustion test facility. 

The structure of the measured spectra differed from the spectra predic- 
ted on the basis of longitudinal duct resonance mode theory using the 
combustor duct system exit temperature to determine the acoustic isen- 
tropic propagation velocity. However, the changes in the spectral 
structure when the g^‘ometry is changed longitudinally caused most inves- 
txg.a'tors to conclude that the spectral structure is related to longitu- 
dinal duct resonance modes. 

Since the spectral structure could not be understood, doubts arose 
about the usefulness of combustion data taken on a liquid fuel open flow 
combustion test facility. The money saved in using this type of facil- 
ity for detailed study of combustion noise would be wasted if the ex- 
perimental results could not be understood. 

The work discussed herein is part of a research program at the 
NASA Lewis Research Center. The purpose of the research program is to 
analyze the pressure spectra measured within an operating ducted gas 
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turbine engine coabuetor test sppemtus. Ibe research prograa has the 
follovlng major goals: 

1. Identify the combustion noise source spectra 

2. Identify the influence of the duct geometry on the spectr'i 

3. Identify Interactions between the combustion process and the acous> f 

tics. j 

This dissertation has the more limited goal of modeling the opera- | 

ting ducted conbustion system and calculating pressure auto-spectra and 
cross-spectra. Consequently, only tvipics related to the propsigation of 

i 

combustion noise in a ducted system are studied. j 

In Chapter II, the failure of the use of the isentropic sound 
propagation velocity to predict the frequencies of longitudinsd duct 
resonance modes in an operating ducted liquid fuel combustion system is 
explained. It is shown that liquid hydrocarbon droplets or solid soot 
particle emissions from a combustor may moke the plane wave propagation 
velocity less than the i.sentropic sp-eed of sound and may attenuate core 
noise. Any variation in plane wave propagation speed may have impor- 
tant consequences in the oesign and selection of acoustical duct lining 
to attenuate engine core noise. This variation in plane wave propaga- 
tion velocity also implies that in order to predict the frequency var- 
iation of resonance caused peaks and dtps in a combustion noise spec- 
trum measured in a duct, one must measure soot and hydrocarbon emis- 
sions and understand the physics and chemistry of the acoustic emission 
interaction. 


In the cases studied, the air mass flow rates are less than 
Kg/ sec, the duct exit temperature is Or’O K, and the flow Mach number 
is much le.ss than unity. The type and amount of emissions depends on 
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such combustor operating conditions as the air and fuel mass flov rate 
and on the type of fuel. It also depends on combustor desiff!) of such 
items as svlrlers and fuel nozzles. Since the paroblem is cc^lext only 
two fairly sinqple types of interactions between emissions and arcustics 
are treated herein. The interaction discussed occurs in a mixture of 
an ideal gas and either a cloud of liquid fuel droplets or a cloud of 
solid soot particles. In both cases, the interaction is between a plane 
wave propagating in the mixture and the cloud. 

Smoke emissions are generated by aircraft engines (Jones (l97T)). 
Consequently, the presence of a cloud of solid soot particles in a liq- 
uid fuel ducted combustor system Is possible. Thus the selection of 
this case as one to be studied seems reasonable. The selection of the 
situation where the duct contains a cloud of liquid fuel droplets is 
more complicated. 

In addition to smoke, aircraft engines also produce gaseous hydro- 
carbon pollution (Jones (1977)). Moreover, Jones (1977) says that un- 
der certain conditions whert' a combination of factors produces large 
dn^ps wd poor mixing of the air with the fuel spray high levels of hy- 
'ii'ocarbons can be im'duced. The hydrocarbon gas pollution may be due 
to droplets which pass through the flame front without oxidizing and 
which eventually eva|K>rate. Consequently , h cloud of liquid fuel drop- 
lets may exist in the ducted combustion system studied if droplets can 
pass throiigh the flame front and if a typical droplet does not evapor- 
ate before It leaves the duct. For the cases studied a particle is in 
the duct less than 0.3 second. The calculation of the evaporation rate 
of a liquid fuel droplet is discussed by Graves and Bahr (1959), Faeth 
(19TT), and Harrje (19TC). Calculations are done using simplified 
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model n. Tlte reaults aro a function of the Initial conditions sueh aa 
droplet dl»u«etcr. Kxaminatlon of nome typical results showed that while 
the ambient temperature (920 K) In JOO to *»00 K above the distillation 
temperature of the liquidy conditions could exist such that a liquid 
droplet islght not evaporate completely in the tiste it is in the duct 
(Uraves and Bahr (1959)). Accordingly, calculations were made consid- 
ering the presence of a cloud of liquid fuel droplets. 

The literature contains a number of theoretical and experimental 
studies of the propagation of a plane wave through a cloud of particles 
which show that particles can cause dispersion wid attenuation. Stud- 
ies considering viscous and thermal interaction but not mass transfer 
were made by Marble (1969), Epstein and Carhart (1953), and Chow 
(196^), Temkin and Dobbins (1966), and Dc-Mns and Temkin (1967). 

Studies which consider water vapor mass transfer in addition to vis- 
v'ous and theimal Interaction were made by Cole and Dobbins (1970), 
Marble «uid Woolen (1970), Davidson (1975), and Marble and Candel (1975). 
Tlie theory of Cole and Dobbins was experimentally conflnned by Cole 
and Dc>bhlns (1971). 

The pliyslcs of liquid hydroearbon droplet attenuation is assumed 
to be identical to the physics of water dmplet attenuation since the 
same physical models ai-e used to describe droplet evaporation (Faeth 
(197TU. The approach followed in caloulatinp liquid hydrocarbon 
vlrx'qvlet att(»miation is similar to that of Marble and Candel (197*’). 

Tlie im'del for the physics and chemistry of carbon particle attenuation 
is h.'uidlod V'.v n mv’di fi cation of the theory given by M n-ble and Candel 
(1Q7‘'). The cons i derat ion of mass trfuu'-fer from the s.oot particle to 
the gas by particle oxidation is new. 
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In the approach used, the particle spacing is assused to be saall 
enout^ in comparison vith the acoustic vave length that the particle 
phase may be represented as a continuous distribution of weak acoustic 
source and dipoles. This is the approach used by Marble and Candel 
(1975) and a similar approach vas previously used by Temkin and Dobbins 
(1966), Marble and Wooten (1970), and Davidson (1975)* The method is 
popular since it shows explicitly the relaxation character of the prob- 
lem and includes the dispersion of sound. 

In Chapter III, a linear analytic procedure for modeling the com- 
bustion duct system acoustics is discussed. The model includes both 
the source region governed by the solution to a non-homogeneous wave 
equation and the propagation region governed by the solution to the 
homogeneous wave equation. All non-linear processes in the system siu:h 
as those which might arise due to interaction of the pressure wave and 
the emissions are neglected in the model. 

A solution to the Sturm-Liouville non-hoiK)geneous one-dimensional 
wave equation, which governs the source region When there is no mean 
uniform flow, can be found by the n«thod of eigenfunctions (Morse and 
Feshback (1973, p. 793), Byron and Fuller (1970)), or by construction 
of the one-dimensional Green's function (Morse and Feshback (1973, 
pp. 122 , 523, 825, 828, and 832), Duff and Naylor (1966)). Applications 
of both solution approaches to aero-acoustics are discussed by Goldstein 
(1976). The eigenfunction solution method was applied to a combustion 
source problem by Mailing (1963). The hon-homogeneous one-dimensional 
wave equation cannot be class! I’ied as being a Sturm-Liouville type if 
flow is included since this introduces complex coefficients into the 
Wave equation. In the analysis developed in Chapter III, the presence 
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of a unifom mean duct flow and the interaction of the flow with com- 
bustion emissiouB is included* and the solution approach is based on 
construction of the one-dimensional Green's function. 

A solution to the homogeneous one-dlMnsi^ial wave equation which 
governs the propagation of plane waves in the duet is formulated as a 
tvD-port transmission line using acoustic transmission matrices. Tte 
vise of this approach to model distributed parasieter systems is dis- 
cussed by Takahashi* Robbins* and Auslander (1972)* and an introduction 
to its use in acoustics is given by Lampton (1973). The equations de- 
veloped by Lan?»ton (1978) are for the no-flow case. Transmission ma- 
trices which include flow have been previously applied to duct acoustic 
problems by Parrott (1973). 

In Chapter IV, the experimental results are presented. In Chap- 
ter V, the ej^erimental results are con^iared with the theory developed 
in Chapter III. 

In sisnmary, this dissertation analyzes dispersion in an operating 
liquid fuel ducted combustion system, gives a procedure for calculating 
pressure auto-spectra and cross-spectra, presents measured auto-spectra 
and cross-spectra, and compares the measured spectra with the calcula- 


ted spectra. 


II. PHYSICAL OONSISQUriOXS 


In this chapter the propagation of a plane vave through a station- 
aiy gas containing a particle cloud of either liquid fuel droplets or 
solid soot particles is discussed. However* the governing equations will 
be derived for the sore general case of non-zero swan velocity since the 
next chapter considers this case. At an appropriate point in this 
chapter the additional aissuoption of zero sietu) velocity will be made. 

A. Governing Equations 

The propagation of a plane wave through a gas containing liquid 
water droplets or solid soot pcu^icles is governed by a wave equation. 
The wave equation used in this chapter was derived by considering a 
plane wave propagating through a droplet or particle cloud in a sta- 
tionary gas. This vave equation was derived by Marble and Candel 
(1975) for a cloud of water droplets. This section reviews the deri- 
vation and interpretation of this wave equation. 

The wave equation derivation is based on the following major as- 
sumptions: 

(1) The particulate mass fraction and volume fraction are small. 

(2) The particulate spacing is smaller than the acoustic wave- 
lengths considered. Consequently, the particle cloud and gas can be 
treated as a continuum. 

(3) The fluctuations of pressure, density, temperature, entropy, 
and velocity are small compared with their mean values. Therefore, 
the squares and cross products of these parameters can be neglected. 
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(k) ']%« fluid velocity is such saaller than the speed of sound. 

(5) The hulk gas can he treated as an ideal gas. 

In deriving the one^diaensional, constant area continuity equa» 
tion» the liquid droplets or soot particles are assuswd to be mss 
sources. Consequently, the density of the gas and the gas velocity are 
related in the continuity equation 

3p/36 + 3(^)/3x « n? (2.1.1) 

where p is the gas density, u is the gas velocity, 6 is tine, x is 
a Cartesian coordinate, and n^ is the local mass production rate due 
to the presence of n particles per unit volume each having a rate of 
mass production of (The symbols used in this chapter are defined 
at the end of the chapter.) 

The momentum equation for the gas is 

3(^)/30 + 3(pi?)/3x + 3p/3x = nF + n^Up (2.1.2) 

where p is the gas pressure. The first term on the ri^t side of 
Eq. (2.1.2) represents a volumetric force due to the n particles per 
unit volxame each exerting a drag force on the gas, F, as the cloud of 
of particles moves through the gas. Using Stokes' drag law, the drag 
force is 

F = 6vru(Up - u) (2.1.3) 

where r is the particle radius, p is the gas dynamic viscosity, and 
Up is the particle velocity. The second term on the right side of 
Eq. (2.1.2) represents the momentum added to the gas locally which is 
equal to the product of the mass generated and the local particle vel- 
ocity u . Substituting Eq. (2.1.1) into vq. (2.1.2) yields 


( 2 . 1.0 
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p3u/»e ♦ pu 3u 3x ♦ 3p/3x ■ nF > - u) 

The ges c<»telnltig the perticles ie esswed to be « aultleoapoaMt 
ideal gas ehareeteriied by the foXloving two equations: 

p - piFt/OW) (2.1.5) 

and 

e ■ c t (2.1.6) 

qp V 

where 5? is the gas c<Mi8tant, t is the tesq^rature, (MH) is the gas 
aolecular weight, e is the interaal energy per unit aass, and is 
the gas specific heat at constant volume. The following entropy repre- 
sentation of the equation of state Is obtsdned by direct integration of 
the relation 

t ds » de - p dp/p^ (2.1.7) 

after sxAstitution of Eqs. (2.1.5) and (2.1.6): 

s - 'So « c^ln[(p/p^)(p^/p)''^] (2.1.8) 

where 

Cp « Cy ♦ a/{W) (2.1.9) 

and 

Y » Cp/Cy (2.1.10) 

where s is the entropy of the gas, Cp is the gas specific heat at 
constant pressure, and y is the specific heat ratio of the gas. 

Fluctuations are assumed to be small so that each quantity is only 
slightly perturbed. Hence, each quantity can be written as the sum of 


n 

an unpertiu1>cd or constant value designated the subscript o and a 
small perturbation designated by the subscript 1: 


P ■ Po ♦ Pi 

(2.1.U) 

U • Uq ♦ u^ 

(2.1.12) 

^ «p.l 

(2.1.13) 

P • Po ♦ Pi 

(2.1. lU) 


(2.1.15) 

♦ “ ♦o " ♦! 

(2.1.16) 

F ■ F + F, 
o 1 

(2.1.17) 


Substituting Eqs. (2.1.11) to (2.1.17) into Eqs. (2.1.1), (2.1.U), 
and (2.1.8) yields the following first-order system of equations: 

DPl/De + 3Uj^/3x = n4^ (2.1.18) 

PqDUj^/DS + 3pj^/3x » nFj^ ♦ "♦l^'*p.o " '^o^ ”*o^'^,l “ '*1^ (2.1.19) 

D(Sj^/Cp)/D0 * D(Pj/yp^)/D 0 - D(Pj^/p^)/De (2.1.20) 

where D/D0 is the substantial derivative, D/D6 ■ 3/36 u 3/3x. 

o 

Equation (2.1.20) is substituted into Eq. (2.1.18) to remove the 
(tensity as a variable. The resulting equation in terms of non-dimen- 
sional parameters is 

D(Pi/yp^^)/D 0 + 3(uj/c^)/3x » n^^/p^ + D(s^/Cp)/De (2.1.21) 


where c is the isentropic speed of sound given by 
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■ Y#t /(lif) • YP«/Pft (2.1.22) 

o o o o 

Equation ( P.1.19) bo vritten In teim of nondiaensional paraMtera 
as 

D(Uj^/Co)/D6 ♦ Cq a(Pi/YPo)/ax - nT^/p^e^ ♦ 

* ( 2 . 1 . 23 ) 

The vave equati«i disctttsed in this chapter is deri>cd by asking 
the following additional assuaptions: 

(1) The average mass production is zero. 

■ 0 (2.1.24) 

(2) The Bean particle and gas velocity are both zero. 

Up^o * % ■ ® (2.1.25) 

Consequently, the one-dlaensicmal, mixture, linearized conservation 
equations for zero mean flow relating perturbations in velocity, pres- 
sure, and density arc as follows: 

^{Pj^/yPq)/^^ + Cq 3(uj^/c^)/3x * + 3(Sj^/Cp)/36 (2.1.26) 

3(Ui/Co)/>6 + Cq 3(Pj^/YPp)/3x ■ n^i/PoCo (2.1.27) 

A wave equation solution to these equations is found next. 

The conservation Eqs. (2.1.26) and (2.1.27) contain unknown source 
terms s^, and In order to solve Eqs. (2.1.26) and (2.1.27). 

first, a small perturbation in pressure is assuwd to cause a saall 
perturbation in each source term. Consequently, transfer functions 
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JlMh • AW)h « and J^(w)/p c relating a swdl change in aass, 
o p o o 

entroi^t and drag force can be defined a& follova: 


[n^i/pj - 

(2.1.28) 


(2.1.29) 


(2.1.30) 


where *[ ] is the Fourier tranafora operator and u is the angular 
frequency. 

The next etep is to take the time and space Fourier transfora of 
Eqs. (2.1.26) and (2.1.27) to obtain on algebraic equation. All per- 
turbations are ass\tned to be proportional to exp |( ik^fl - iwd )j . Conse- 
quently, the Fourier transfom of a function g of the independent 
variables 6 and x produces another function denoted 0 of inde- 
pendent variables u and k^Q defined by 


G((i),k G) 
o 



-Ik tlx 

g(0,x)e 


dx de - ♦[g(e,x)] 


(2.3.31) 


The inverse Fourier transform of a function G of the independent w- 
tables u) and k^Q produce another function* g* of the independent 
variables 9 and x defined by 


g(®.x) 




ik (lx 
fl)e ° 


d(k n)dbi 
o 


■ ♦ 


-1 


[o(u»,k^n)] 


(2.1.32) 
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To find tlm Ftouricr trunc rom uf n time derivnti^, the tiae derivative 
of Eq. (2.1.32) ie taken 


ag(e,x)/ae - 


(2.1.33) 



Since the Fourier transform operation is unique, taking the Fburler 
transform of Eq. (2.1.33) yields 

•[3g(e.x)/ae] ■ (-i«)#[g(e,x)] ( 2 . 1 . 34) 

The corresponding Fourier transform of a space derivative is 

♦ [9g(8.x)/3x] ■ (ik^0)s[g(e,x)3 (2.1.35) 

Taking the time and space Fourier transfom of Eqs. (2.1.26) and 
(2.1.27) and substituting Eqs. (2.1.28) to (2.1.30) yields the follow- 
ing equation in matrix form 



(2.1.36) 

A non-trivial solution exists providing the determinant of the coeffi- 
cient matrix vanishes. Solving the determinant yields a second-order 
polynomial equation for (ik^fl) 
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( 2 . 1 . 37 ) 


The solution to Eq. 


{?. 1.37) is given by 



An approximate solution to the wave number equation (Eq. (2.1.37)) is 
given by Marble and Candel (1975). This solution may be obtained from 
the exact solution (Eq. (2.1.36)) by assuming that 


and that 



(2.1.U0) 


( 2 . 1 . 1 *!) 


With theae assumptions and after son^ algebraic manipulation the two 
approximate wave number solutions to Eq. (P.1.37) are 
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(2.1. k2) 


(2.1. J»3) 


The wave Qximber for vaves moving to the right used by Marble and Candel 
( 197 ^) is given by Eq. (2.1.142) in terms of the notation used herein. 

By definition, the sound propagation velocity is related to the 
wave number eqi»rtion by 


Also by definition, the acoustic energy attenuation coefficient in 
Nepers per meter is 

a = -2 Re(ikQJ2) (2.1. U 5 ) 

Nepers per meter is converted to decibels per meter by multiplying by 
10 log e = U. 3 I*: 

a = -8.68 Re(ik 0), dB/m (2.1.U6) 

o 

In this section the fundamental governing equations were discussed. 
Furthermore, for the case where the average mass production , mean parti- 
cle velocity, and mean gas velocity ;vre zero, equations for calculating 
the sound propagation speed and the attenuation of a plane pressure wave 
propagating in a cloud of peurticles were determined. These equations 
depend on transfer functions relating a small change in mass production, 
entropy production, and drag force to a small change in pressure. The 
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derivation of these transfer functions for the case of a plane vave 
propagating in a cloud of liquid fuel drops and a clor^ of solid soot 
particles Is discussed next. 

B. Response of the Fuel Ihroplet Cloud 

Fuel droplet ealsslons ttcm coabustors in the c<»bu8tion duet sys- 
tem can occur when the fuel droplet spray consists of large-dlaaeter 
drops and the air Injected around the fuel notsle through svlrlers 
mi <s8 poorly with the fuel spray. This combination of large fuel drops 
and poor mixing produces high levels of hydrocarbon emission (Jones 
(1977)). The larger droplets which pass throu^ the high temperature 
reaction sone are assumed to evaporate iu a relatively cool gas mixture 
having a low oxygen concentration. Accordingly* the process is approx- 
imated by considering evaporation without combustion (Faeth (1977)). 

The processes and assumptions typically employed in aui analysis of 
fuel drop evaporation are discussed by Graves and Bahr (1950), Faeth 
(IQ"^?), and Harrje (1972). However, rather than using a complex model 
to obiain expressions for v#(w)/p^, J(w)/Cp, and .^(«)/ppC^, the 
approximate model used by Marble and Candel (1975) for water drop evap- 
oration is adapted for use herein by making appropriate assumptions. 
'PiUs model is discussed in Appendix A. Tlie resulting expressions pro- 
vide: (l) a guide to the order of magnitude of the transfer functions 

.4f(u))/p , J(u)/c , and c , and (2) Information on the fre- 

o Y o o 

quency dependence of the soimd propagation speed and the attenuation, 
because the resulting expressions are identical to those which cotUd be 
obtained from expresaiona derived by Marble and Candel (1975), they are 
shown in Appendix A. The numerical values of the factors in these ex- 
pressions are different because the physical properties were changed 
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to thoie for o fuol drop mther than thoit for m wotor drop. Tlio ro« 
suits of eolculotions using this theory ere discussed in section D. 

C. Response of the Soot Perticle Cloud 

Soot produced hy ges turbine cosibustors oxidises es it flows 
through the duct systea (Linten end Heyvood (1971) » Rorgren (1971) • end 
Oouldin (1973)). In order to eccuretely celculete the response of e 
cloud of soot perticles to en ecoustic pressure perturbetion. the aess 
distribution, the chesdcel coaposition , end the distribution of perticle 
shapes and sises should be known. It is also necessary to know the rate 
constants for the heterogenous oxidation process and the heat transfer 
process. However, at this tine none of this information is known ac- 
curately. Conseqxwntly, the following model includes the basic physics 
and chemistry in only an approximate manner. An experiTOntal study of 
the oxidation of soot by Wright (1975) and a mathematical model for the 
gasification of coal char discussed by Siisaons and Lewis (1977) were 
useful in formulating the approach used in this work. The model is 
based on the assumption that the soot particles can be treated as being 
made of porous carbon, having a spherical shape, and having a uniform 
temperature. 

The transfer functions »4f(u)/p^, J(u)/Cp, and J^(u)/p^c^ are 
derived from equations that describe the particulate phase of the one- 
dimensiona] flow of the gas-soot particle mixture with transfer of 
mass, momentum, and energy. When specifying the functional dependence 
of and ,l(w)/c„, the following selection is made. As pre- 

viously diiicussed* for the force acting on the gas due to a particle, 
Stokes’ viscous drag law is used. Thus, the total force on the gas is 

F * -n6i^r y(u - u ) 
s s 


( 2 . 3 . 1 ) 
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For the heat transfer froa the particlee to the Hewton's lav of 
conirection is itsed. Thus* the total heat transfer frcns n soot parti* 
cles Into a unit vol\a»e of gas having thersMl conductivity k is 




?^U(t..tj 


(2.3.2) 


where the heat-treinsfer Nusselt number NUj^ equals 2 if the sdr is 
8ta®iwrit. Tbe value 2 which is generally assvmed, is used herein since 
the difference between the particle velocity and the gas velocity is 
assumed anall. 

The viscous drag transfer function is discussed first. A soot 
particle of mass m and radius r is moving with a velocity u . 

S @ S 9 JL 

in the gas. The gas has a velocity u^ due to the wave motion. The 
particle experiences a Stokes* drag force of 6srgU(u^ “ % 1^’ 
resulting particle equation of motion is 


"s '■'l ■ - “s,l’ 


(2.3.3) 


The gas velocity is related to a pressure perturbation for small per- 
turbations by 




(2.3.1*) 


The time Fourier transform of Eq. (2.3.3) yield.*? the following expres- 
sion for the response of the particle velocity to a pressure perturba- 


tion: 


■ rT-ferr *[^] 


(. 5 . 3 . 5 ) 
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^ere 


T ■ 


S 


(2.3.6) 


Physicallyt the tiae ooastent t is the tiae interval in vhich the 
particle velocity is reduced a factor of l/e froa its initial 

value due to tl» Stokes' drag force. Substituting Eqs. (2.3.U) and 
( 2 . 3 . 3 ) into Eq. ( 2 . 3 . 3 ) and using Eq. (2.1.30) yields the folloving 
expression for the viscous drag force perturbation in sresponse to a 
pressure perturbation: 


n6vr u c (-iw)T 
s 0 

A 

^1 

PqC^ 1 + (-i«)t 




(2.3.7) 


Consequently, the viscous drag force transfer function Is given by 




1C 

s 

(-iwh _ 

T 

_1 ♦ (-iw)x_ 


( 2 . 3 . 8 ) 


Where 


•c 


s 


nm 

s 



( 2 . 3 . 9 ) 


The heat generated by the oxidation of a soot particle Is assumed 
to be transferred directly to the particle. Thus, the entropy source 
perturbation term may be calculated considering only the heat transfer 
between the n particles and the bulk gas using Eq. (2.3.2) as follows: 


, “'1 

^ P 

00 d 6 








) 


(2.3.10) 
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The differential equation for the entrojor source perturl>ation ean then 
be vritten 


A 

dO 



( 2 . 3 . 11 ) 


yhere 


T 


8 


m c 



(2.3.12) 


Physically, the tine constant Is the time interval in which the ^ 
entropy is changed by a factor of 1/e from its initial value. 

Perturbations in bulk gas tenperat\ire and pressuire are related to 
perturbations in soot particle tenperature by combining the tine deri- 
vative of the bulk gas energy equation (Eq. (2.1.20) with u^ » 0.0) 
with the perfect gas law and Eq. (2.3.11) to show 


_d. 

de 






(2.3.13) 


The slowest process in the burning of a soot particle controls the re- 
action rate. Depending on the temperature and the particle geometry, 
this may either be the diffusion of oxygen to the particle siorface or 
the chemical reaction rate at the surface. For this model, the chemi- 
cal reaction rate is assumed to be slower than the diffusion rate 
(Linden and Heywood (1971)). Conseqwntly, it controls the oxidation 
process. The soot is, for simplicity, assumed to be carbon. However, 
evidence exists that the soot contains a few percent hydrogen which, 
in practice, could be important. To model the chemical reaction of 
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oxjrsen vith the soot» an adsorption/desorption aechanisa la uaed in- 
volving the follovlng sequence of reactions: 


c ♦ Og - 

X 

“ (“a). 

(“2).- 

'a, 

fc. C ♦ O2 



K- 

c + (CO2' 


2 C 0 


where ^^ 2 } ^ surface oxide (i.e.» an O 2 loolecule adsorted on the 

carbon surfaces (Wright (1975) » Simmons and Lewis (1977)). 

This process may be followed by a gasification reaction which pro- 
duces carbon dioxide (Wright (1975)). Only the initial production of 
carbon monoxide is considered herein. 

The following concentration eqmtions can be written for the se- 
quence of chemical reactions considered: 


* ‘ ''iN ■ "2 [(“ 2 ).] • S [(“s).] 

- = -■'2 [(“s).] ’ ■'iN 

The concentration of ( 0 O 2 ) is assumed constant because a steady-state 

' ^a 

condition is achieved where as much is created eis is destroyed. Conse- 
quently, the time derivative of is zero. Thus, the left side 

ft 

of Eq. ( 2 . 3.1^*) can be set to zero, and the following expression for 
the concentration of ^^© 2 ^ is obtained: 



1 K, c 

K)J • ^ 




(2.3.16) 


The cerbon consunqption rate» vhich equals the net rate at vhleh aass is 
added to the gas by a single particle per unit particle surface area, is 
then given by 








Kg ♦ Kg 0t 


( 2 . 3 . 17 ) 


* Vo, 


The following Arrhenius model chsnlcal reaction rates (from Simons 
and Lewis (1977)) were determined empirically to fit a wide range of 


data: 


= 9?00 e“19,700/t 


( 2 . 3 . 18 ) 


K,/^t » 2xl0”^ ^-15»700/t jjg/jj_ggg 


(2.3.19) 


K, . 5.5 10^ e-3e-7<>0/‘ 


(2.3.20) 


The rate of mass addition to the bulk gas per unit volume for soot par- 
ticles having surface area is then 


= n^_ + 





( 2 . 3 . 21 ) 


For a small perturbation in surface temperature and oxygen particle 
pressure, the resulting change in mass addition can be written as 


= na K 
s a 
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Ibus, the otse perturhetion term is given hgr 


n#, ■ (na 


K -n )( A Aq S.l 


(2.3.23) 


This equation ie rewritten as 






°2 PO, 


2.0 


(2. 3.2‘*) 


where 


T B i— . 


A 


I 


a 

s 



Y 



and 


C 


t 

o 

K 

a 



(2.3.25) 


(2.3.26) 


(2.3.27) 


( 2 . 3 . 28 ) 


The external sxirface area per unit mass of a spherical carbon particle 


of diaxneter d is 

5 





m 


s 



1* /d\^ 2000 ^ 

3 *^2; 


1000 d 


hg 


(2.3.29) 


Typical internal surface area per unit mass of coal char as measured 

2 

by molecular adsorption are 100,000 to 1+00, 0(X) m /kg for 100 w particles 
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which ia four orders of magnitude greater than the external surface au?ea 
per unit aasa (Siamons and Levia (1977) )• Conaeq\i»ntly, the model cal- 
culations made Iwrein use internal surface areas per unit mass approxi- 
mately four orders of magnitude greater thu the external surface areas 
per unit mass. 

The oxidation reaction is 


C ♦ I Og CO 


Ihe heat generated is given hy 




(2.3.30) 


^e heat balance equation for a soot particle is then 


(^) '“..1 ■ '-. 1 > - Vl 


When all n particles are considered* Eq. (2. 3-31) becomes 


s / s.i - 


(2.3.32) 


where 


(2.3.33) 


The perturbation in the partial pressure of oxygen is related to the 
Bass source variation and the total pressure perturbation using the 
continuity equations. The oxygen continuity equation is 
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36 





-n6 


1 12 


(2.3.3>») 


for when a aasti of carbon, 6^, in added to the gas ttie oxygen concen- 
tration is reduced by l6 6/12. Tbe bulJt gas c<»tia\dty eqnstiM Is 


i-p 4p 
36 *^1 *^0 3x 


n6. 


(2.3.35) 


The velocity gradient is elisdnated from Eqs. (2.3<3^) snd (2.3.35) ss 
follows : 


Thus, 


3x p^ 


Po 


-M 





(2.3.36) 



(2.3.37) 


However, the oxygen and the bulk gas temperature are identical, and 
the gas properties are assumed to be identical; that is, both behave 
as perfect gases. Consequently, Eq. (2.3.37) can be written in terms 
of pressure perturbations as follows: 



Equation 
(2.3.32), and 


(2.3.2*<) and the Fourier transforms of Eqs. 
(2.3.38) can be put in the following matrix 


(2.3.13), 

form: 



( 2 . 3 . 39 ) 



3k 


The solution Is 

V ■ (2.3.J»0) 

The sMiss source transfer functicm is found hjr using the solution 
tor Ij/pQ tram Eqs. (2.3*3S)) and (2.1.28). The entropy 

source transfer fimction is found using the quantities • [t -./t^l and 

8 gX O 

♦ [t_ ,/t^] which were obtained from Eq. (2.3.39)* the Amrier trans- 

*x o 

fora of Eqs. (2.3.U) and (2.1.29). 

In this section the response of a soot particle cloud to a pres- 
sure perturbation was investigated. As a result equations were derived 
for calculating the transfer functions relating a small change in mass 
production, entropy production, and drag force to a small change in 
pressure. The next section discusses typical dispersion and attenua- 
tion curves calculated for sound propagrtion in a liquid fuel droplet 
cloud and in a soot particle cloud. 

D. Sasq>le Calculations 

A typical dispersion curve and an attenuation curve calculated for 
vaporizing fuel droplets by using the arbitrary parameters in Table I 
are presented in Figure 2. The propagation velocity varies from 506 
m/sec at low frequencies to the isentropic speed of sound (610 m/sec) at 
high frequencies. The attenuation is greater than 3 dB/m above UOO Hz. 

The relative inqxjrtance of the transfer functions Jf[u)/p . 

o 

and e in calculating attenuation and dispersion 

pop 

using Eqs. (2.1.42), (2. 1.1*4), and (2.1.46) can be determined trtm Fig- 
ure 3. The viscous drag transfer function, c , and the mass 

o o 

source transfer function, u^(“)/Pq. shown in Figure 3 have negligible 
effects at combustion noise frequencies which are below 500 Hz. Conse- 
quently , the entropy source transfer function determines the attenuation 
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TABLE I. - P*r«seters^ used to calculate dispersion and attenuation 
for the case of fuel droplets vaporizing 



. . . . 922.0 

V 

.... 0.05 

K, W/m-k . . . 

. . 5.38/100"® 

«VL- 

. . k-iexio"^ 

Pp. kg/m^ . . . 

. . . . 0.378 

nBj^, kg/m 

. . . 0.0058 

Cp, J/kg-K . . 

. . . . 1100.0 

c^, m/sec 

• • • • ^XOs 0 

Y 


n, Hy^/T^Cp . . . . 


kg/m^ . . . 

. . . . 81 * 5.0 

(MM)^:W)g. . . 

.... 3.9^ 

J/kg-K . . 

1700 

n 

. . 1 . 6 i*xl 0 ^^ 

“ 

. . 1100 >‘ 10 '^ 

kg 

. 3 . 5 l*>‘ 10 "^^ 

(MW)g 

. . . . 28.97 



, . . 1.53*10"® 

<»“)l 

. . . . lll ».0 

T, sec 

, . . 5.13*10"^ 

W, kg/m-sec . . 

. . 3.66x10"^ 

Tp, sec 

, . . 1 . 6 Uxio"^ 

2 

m /sec . . 

. . I*. 5 l*xl 0 "^ 

^TL* 

. . . 1. 18*10"^ 


^Defined in Appendix A. 




FREQUENCY, Hz 

Figure 2. - Calculated acoustic attenuation and dis- 
persion in region containing valorizing fuel droplets. 
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Figure 3. - CiteuWed viscous drag force, mass 
source and entropy source transfer functions in a 
region containing vaporizing fuel droplets. 



and dispersion at lov frequencies. Specifically* the viscous drag force 
transfer function can be neglected belov 1200 Hz. Above 1200 Hz the 
entropy transfer function becanes less than the viscous drag force 
transfer function and this accounts for the increase in the attenuation 
shown in Figxure 2 above 1200 Hz. 

A typical dispersion curve and an attenuation curve calculated for 
oxidizing soot particles by using the arbitrary paraaeters in Table II 
are shown in Figure for soot particle mass inactions of 0.01^ (exam- 
pie 1, Table II) euid O.OOll* kg/m (example 2, Table II) and an oxygen 
partial pressure of 0.05. The sound pi-opagation speed varies from i»l«0 
m/sec for a maes fraction of O.OOlU and from 1*20 m/sec for a mass frac- 
tion of 0.015 at low frequencies to the isentropic speed of sound (6l0 
m/sec) at high frequencies. The attenuation for a mass fraction of 
0.015 is greater than 3 dB/m above 1*00 Hz. However, it drops to less 
them 1 dB/m for a mass fraction of O.OOll*. It was also found that for 
these cases the viscous drag force transfer function and the mass 
source transfer function cotild be neglected at low frequencies. 

In this chapter the ;ropagation of a plane wave through a station- 
ary gas containing a particle cloud of either liquid fuel droplets or 
solid soot particles was studied. The fundamental governing equations 
were discussed. In addition, sample ctLlculations were made producing 
typical dispersion and attenuation curves for both cases. Sample trans- 
fer function calculations for the liquid fuel droplet cloud case were 
also presented. The failure of the use of the isentropic sound propa- 
gation speed to predict the frequencies of longitudineQ, duct resonance 
modes is explained as being due to the presence of clouds of particles 
which cause the sound to propagate at a lower velocicy at combustion 
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TABLE II. - Paraaeters used to calculate dispersion and attenuations 
curves for the case of soot particles oxidizing 


P JC • • • • 


K, V/b-K . . 

. . . 5.38X10"^ 

3 

p, kg/m'^. . . 

0.378 

Cp, J/kg-K . 

. . . . 1100.0 

y 

l.J* 

3 

Pg. kg/m . . 

. . . . 1880.0 

3 

c„, J/kg-m . 

3 

. . . . 3200.0 

X* a IB • • » • 

S 

. . . l.OxlO"^ 

(MM)g .... 

28.97 

y, kg/m-sec . 

. . . 3 . 66 x 10 "^ 

'°2 

0.05 


Q. J/kg . . . 

. . . .-3.637x10^ 

0 

Aj, m /kg . . 


c .*n/sec . . 


0 

€ % Q/T c • • 


0 p 

K , kg/H-sec 

Ct 

. . . 8.08x10“^° 

T dK 

0 a 

^ k" IT * • 



m , kg . . . 

s 

. . . 7.87x10“^^ 

T, sec . • • 


T a sec a • a 

.... 3.50x10"^ 

S' 

T a sec a . . 


a 


Example 


nm , kg/m 


"s’ 


7.36>‘10 

5.8x10"" 


i.5>'10 


6.73>‘10 

— 

5.3x10 

l.l*xio“^ 







Figure 4. - Attenuation and (Sspersion for oxicSzing 
soot particles. 



noise frequencies. In the next chapter the propagation of a plane wave 
through a noving gas containing a soot particle cloud is discussed. 
Moreover* a Biodel for the operating ducted c^bustor syst«s vill be 
described lAich can be used to calculate in-duct pressxire auto-spectra 
and cross-spectra. 

Symbols Used in This Chapter 
j/ matrix 

2 

Ag soot particle exterior area per tinit mass, m /kg 

2 

Aj soot particle interior area per unit mass* m /kg 

2 

a^ soot particle surface area* m 

c soot mass* kg 

c((«>) sound propagation speed* m/sec 

isentropic speed of sound, m/sec 
Cp gas specific heat at consteuit pressure, J/kg-K 

Cg soot particle specific heat* J/kg-K 

gas specific heat at constamt volume, J/kg-K 
d diamter* m 

e interned energy, J/kg 

F particle force exerted on gas 

^(‘»))/p^c^ viscous drag transfer function 
f frequency, Hz 

2 

H heat-transfer coefficient, W/m -K 

1 (- 1 )"'^ 

K effective chemiceil reaction rate, kg-N-sec 

K, chemical reaction rate of the process 


propagation wave nximber* u/c , m” 

o 


(MW) 

m 

n 



T 

J (w)/Cp 

s 

t 

• 4 * 

U 

u 

X 

-¥ 

Y 

Y 


a 

e 

y 

t 

c 


mass sourc#? transfer function 

molecular weight 

soot particle mass* kg 

heat-transfer Ni^selt number, Hd/r 

number of particles per milt volume 

number of moles 

pressure, N/m^ 

2 

partial pressure of oxygen, N/m 

heat transferred to gas from particles hy convection 
heat generated by chemical reaction per unit mass of fuel, 
J/kg 

gas constant 

particle i*adius, m 

entropy source transfer function 

entropy of gas, J/kg-K 

teagjerature , K 

system forcing function 

velocity, m/sec 

Cartesian coordinate 

system state vector 

relative concentration, partial pressure ratio or density 
ratio 

acoustic attenuation coefficient, dB/m 

see Eq. (2.1.39) 

specific heat ratio of gas 

chemical heat parameter, Q /c t 

c p o 

chemical reaction rate parameter, (t /K ){dK /dt) 

o a a 



I 

I 

I 


i 

1 




k3 


% 


K 

K 

y 

p 


T 

T 

T 


a 

s 


♦I ] 

♦ 
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Superscripts 

T) 

n 

( )* 

( )■ 
a 

C 

H 

°2 

0 

P 

s 

1 


tine, sec 

gas theraal conductivity, W/i^K 
soot particle SMSS fraction, nm^/p^ 
gas dynaaic viscosity, kg/n-sec 
gas density, kg/a 

Stokes* relaxation time, B/6vr.)i, sec 

adsorption relaxation tine, 

soot particle thermal relaxation tisa, 

2 

(mgCp)/lmr (HUjj/dg)ic, sec 
Fourier transfom operator 
mass source rate 

dispersion-attenuation propeigation wave nuatber factor 

angular frequency, 2irf, radians/sec 

and subscripts; 

vector qusmtity 

instantaneous quantity 

downstream propagating 

upstream propagating 

adsorbed 

carbon 

heat transfer 
molecular oxygen 
reference state 
particle 
soot 

perturbed quantity 


in bulk gas far from soot particle 


III. ACOUSTIC MODEL OF THE COMBUSTION DUCT SYSTEM 


The effect of attenuation and dispersion on duet spectra measured 
in a flowing confined gas containing soot particles is discussed in this 
chapter. The analysis presented in the last chapter Investigated atten- 
uation and dispersion in a stationary* infinite bulk gas containing 
particles. Consequently, the previous analysis must be modified to 
take into account the presence of mean flow. 

In the last chapter it was determined that the viscous drag trans- 
fer function .^((o)/p c , and the mass source transfer function, 

o o 

v^(u))/p^, have negligible effects at combustion noise frequencies on 

attenuation and dispersion. Tnus, the entropy source transfer function, 

«^(u)/c , determines the attenuation and dispersion at low frequencies. 

P 

However, the entropy source transfer function, defined in 

Chapter I depends on terms related to the mass source. 

In this chapter, it is assumed at the start that the viscous drag 

transfer function, ^(<«))/p c , and the mass transfer function, 

o o 

uf{ii))/p^, have negligible effect at combustion noise frequencies on at- 
tenuation and dispersion and that only the entropy source transfer 
function, </(u)/Cp, is important. 

In addition, the entropy source transfer function, J(o))/Cp, de- 
fined in this chapter will be determined by assuming that the conse- 
quences of soot particle surface oxidation can be ideEdized as produr 

» 

cing a constant time-independent soot particle temperature. Conse- 
quently, the model presented in this chapter Includes only the effect of 
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heat transfer between the soot particles and the flowing gas. This as- 
suaption yields a set of equations depentent on only a single tine con- 
stant. As a consequence, the equations can be forstulated as a distri- 
buted system parameter identification problem with a single unknown 
time constant. The model described in Chapter I is more coiqilex and 
depends on more parameters. 

An extensive review of the distributed system parameter identifica- 
tion problem is given by Polis and Goodsen (1976). Rather than using 
any of the methods described by Polis and Goodsen to evaluate the un- 
known time constant, it will be evaluated as follows. First, the solu- 
tion to the problem discussed herein is approximated by the solution to 
an acoustic wave equation. Then, the acoustic wave equation and its 
solution are used to model the propsigation in the duct system emd to 
ccLlculate auto-spectra and cross-spectra. Finally, by conq>aring meas- 
ured spectra ai:d calculated spectra the unknown time constant may be 
determined. 

A. Analysis 

1. Model Wave Equation 

The mass transfer and the body force due to viscous drag are ne- 
glected. Consequently, Eqs. (2.1.18) and (2.1.19) yield in non- 
dimensional form the following one-dimensional continuity and momentm 
equations relating small perturbations: 

D(p /p )/D6 + c 3(u /c )/3x = 0 (3.1.1) 

1 o o 1 o 

D(Ui/Co)/De + c^3(p^/yp^)/3x « 0 (3.1.2) 

The energy equation derived by considering only the heat transfer 
between the particles and the bulk gas is 
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tpDs/DO - hQ (3.1.3) 

vhere the heat transfer rate per soot particle froa the i»rticie into a 
unit vol\ae of gas is 

)ie(tg - t.) (3.1. >») 

This equation is similar to Eq. (2.3.10). However, the total derivap* 
tive appears in place of the time derivative since tlM gas has a SKsn 
velocity. The total derivative is given by 

D/D 0 * a/ae t u a/ax 

0 

where the first term on the right side is the temporal member and the 
second term is the convective meniber in the x direction. When the mean 
flow is zero, as it was assumed to be in the previous chapter, the con- 
vective member is zero. However, even with flow, when the entropy is 
constant the convective term is frequently neglected if the time to con- 
vect a disturbance through the volume of interest is large coag>ared with 
the time required for a quantity to change from a minimum to a maximum. 
The entropy is not constant in this problem since the heat transfer be- 
tween the soot particles and gas is taken into account. Therefore, the 
effect on the attenuation and dispersion of the convective term and the 
heat transfer term may be of the same order of magnitude. Consequently, 
the convective tenn is included in the governing equations. 

In deriving the small perturbation form of the energy equation the 
soot particle temperature is assumed to be time-independent 


t » t 
s s,o 


(3.1.5) 


Consequently, for small perturbations the ener^ equation (Eq. (3.1.3)) 


yields 
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D(s,/c_)/De ■ ,/t (3.1.6) 

Thus, the system model depends on a single time constant, 

time constant is the ratio of the heat transfer time constant to the 

soot particle mass fraction. 

These equations are sii^lified further using small pertui^ation 
thermodynamic relations. !I%e idead gas entropy equation is, for small 
perturbations using Eq. (2.1.8), 



(3.1.7) 


Also, using Eq. (2.1.5), for small perturbatiom the gas equation of 
state is 



( 3 . 1 . 8 ) 


Equations (3.1.7) and (3.1.8) are substituted into Eq. (3.1.6) to deter- 
mine the response of the gas tenqjerature to a pressure perturbation. 
Thus, 


M. ^ (II . . J-fti -2l 

D [cj » oj De[ YP„ 



(3.1.9) 


From Eq. (3.1.9) the desired relation is 



( 3 . 1 . 10 ) 


To remove the density perturbation, Eq. (3.1.7) is substituted into 
Eq. (3.1.1). Thus 




kd 



Soiring Eq. (3.1.11) for D(«j^/Cp)/W yields 


D6 



(3.1.11) 


(3.1.12) 


Equs^ioxis ( 3 * 1 . 2 )^ (3«1*6)( (3.1.10)^ snd (3*1.12) sre fozids^ 
mental model equations, niese equations are next vritten in state var> 
iable form using acoustic pressure* particle velocity, entropy, and gas 
tenqperature as the state variables. Then a solution to the model dif- 
ferential equation is derived. 

2. State Variable Formulation 

The system differential equation based on the small perturbation 
approximation is discussed next. Equations (3.1.2), (3.1.6), (3.1.10), 
and (3.1.12) are written in state variable form as 




(3.1.13) 


where 


1 - Y 
0 
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0 

1 

0 
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0 
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1 

-1 


1 

0 
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(3.1.1U) 
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C^/u 
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( 3 . 1 . 15 ) 
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K /t 
S 8| 


(3.1.16) 


’l/l'Po 1 


“/'o 


’IS 


( 3 . 1 . 17 ) 


The state variable solution to Eq- (3. 1*13) is asswed to have the 


4 


(ik fl),x-i(.'e 
o t 


(3.1 .18) 


Taking the Fourier transform in the time variable and the space varia- 
ble of Eq. (3.1.13) produces a set of four homogeneous algebraic equa- 
tions in four unknowns expressed by the following matrix equation 


tl 

p 
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[a(-1w) + B(ik^n)u^ ♦ c]<![( -i«), (ik^Q)] (3.1.19) 


A non-trivial solution to Eq. (3.1.19) exists if and only if the follow- 
ing determinant of the coefficient matrix vanishes: 


A(1 - y) 


A.- 

s 




(3.1.20) 


where 


c (iic Q) 
o o 


IC /t 
s s 


A ~ (-iw) + (iic Q)u 
o o 


Solving the determineint equation yields a third-order polynomial wave 
number equation for (iK^fi), 




and a solution A = 0 which corresponds to a wavenumber 


^ vO 103 
{IK = — 

o a u 

o 


The roots of the third-order polynomial equation can be found numeri- 
cally to great accuracy using an iteration method due to Muller (Conte 
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and de Boor (1972)). For the data studied herein the Nach nunber is 
low and a good approxinate solution at all frequencies is 




9 ^ i<tf 


T U 
8 O 


i(«/c ) 

■ mr 

s 


(iK 0)“ = 

o 


i(w/Cg) 

1 - M_ 


»rtiere 


(3.1. 2U) 


(3.1.25) 


(3.1.26) 



tind 



c 

s 


(3.1.27) 


(3.1.28) 


The state variable solution to Kq. (1.1.13) is expressed in terns 
of four one-dimensional normal inodot? with wave nun»bers given by Eqs, 
(3.1.23) to (3.1.26). Equationr. (3. 1.25) and (3.1.26) indicate that two 
laoden may be interpreted as acoustic modes corresponding to right and 
left traveling acoustic waves. The remixlning inodes represent waves 
traveling with the flow. The wave defii:od by Eq. (3.1.23) is atten- 
uated. However* this is not true lor the wave defined by Eq. (3.1.21*) 
since assumed to be small. Thus* the attenuation of the wave 
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traveling with the flow calculated by substituting Eq. (e.l.2U) into 
Eq. (2. 1.1*6) is very large emd this wave is highly damped for the cases 
considered herein. Consequently, this mode quickly decays in the cohh 
bustion duct system. The wave with wave number given by Eq. (3.1.23) 
is assuoed to have negligible effect on the system. Thus, the pressure 
and particle velocity and the structure of the pressure auto-spectra 
and cross-spectra are assumed to be determined only by the acoustic 
inodes. 

The state variable differential equation (Eq. (3.1.13)) depends 
upon four dependent variables: pressure, particle velocity, entropy, 

and bulk gas temperature. In the next section the non-acoustic modes 
are removed from this state veuriable differential equation. This pro- 
duces an acoustic state veu"iable system differential equation which has 
two independent variables; pressure and particle velocity. However, 
the new acoustic wave number solutions are identical to the previous 
ones. The acoustic state variable system differential equation is 
further simplified by use of a velocity potential function to create a 
single one- dimensional differential equation. This is the equation 
used to model the system. 

3. Acoustic State Variable Formulation 

The set of differential equations which has only the acoustic mode 
is obtained by neglecting the convection term, u^9s/9x in Eq. (3.1.6). 
As a consequence Kq. (3.1.9) becomes 




1- 

(1 - y) Pi ^ 

K t - 
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(3.1.29) 
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l*Yom ihl'A oquailon ihr mni^onno of the bulk p.ac temperature to a pres* 
sure perturbation Is 


/a ^ !3\ . . 4 Il^\ !i 

r "-.0 V Y ho 


(3.1.30) 


Taking the Fourier transform in the time domain in Eqs. (3.1.29) and 
(3.1.30) yields 


, 1~ | _ (-iti>)(l - Y 





and 


s 

c 




YP, 


Where 


iiid 

C 


-0^7. ll. 


1 ^ (-i.) 


The time Fourier transform of Eq. (3*l.l) is 


(-iu) + u 


b 3x 


L o. 



•\ r 1 

+ 0 I = 0 

o 3x YP^ 


i-'ubstitutinfi Eq. (3.1.3.?) into the time Fourier transform of 
Eq. ( 3.1.1;'’) yields 


C • 4. , -J 

, .rcu)] 


r-'il 
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IP 
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J 0^ 
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= 0 


(3.1.31) 


(3.1,32) 


(3.1.33) 


(3.1.314) 


( 3 . 1 . 35 ) 
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A plane wave solution to Eqs. (3.1.3^) and (3.1.35) can be obtained In 
terms of a velocity potential, where 

(3.1.36) 

and 


Pi ■ "o (m * “o >i)* <3.1.37) 

Substituting the tine Fourier transform of Eqs. (3.1.36) and (3.1.37) 
into Eq. (3.1.35) yields a new governing partial differential equation 
called the non-adiabatic velocity potential wave equation 




r 

(-‘Xs’ * ij 


■‘Xs> * «3 if*'*' - 7? 

dx 




(3.1.38) 


Where 


c 

s 



(3.1.39) 


and 


k ^ 



( 3 . 1 . 1 < 0 ) 


The velocity potential solution is assumed to be proportional to 
exp((ik^fi)x - ioO)). Substituting this solution into Eq. (3.1.38) 
yields a wave number equation which has the acoustic wave numbers given 
by Eqs. (3.1.25) and (^.1.26) as a solution. Consequently , the velocity 
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potential wave equation solution consists of an acoustic vave traveling 
upstream and downstream expressed as 

r (ik n)^ (ik nTx"! 

♦ - ae ° ♦ be ° Je (3.1.1*!) 

In the following sections the use of the Fourier transform operator no- 
tation is discontinued and all variables such as p, and u are 
Fourier transformed variables. The following sections apply the acous- 
tic state vector formulation and the velocity potential definition to 
the study of sound propagation in the ducted combustion system. 

B. Application to Ducted Combustion System 

The analysis developed in the last section is now used to study the 
spectral structure of pressure measurements made in the ducted combus- 
tion system shown in Figure The ducted combustion system sho%m in 
Figure 'j> consists of: (l) a source region inside the combustor can, 

(?) a non-source region inside a spool piece and a long duct, (3) an 
area expansion and contraction on either side of the spool piece, and 
(I 4 ) a downstream boundary at the exit of the long duct and an upstream 
undary at the combustor inlet. 

The solution for the velocity potential in a non-source region is 
given by Eq. (3.1.!4l). The velocity potential wave equation 
(Eq. (3.1.38)) is assumed to apply in the source region with the addi- 
tion of a source term G(u,x) on the right-hand side. Using the com- 
bustor inlet and exit impedance as a boundary condition, a unique Green's 
function solution for the velocity potential is foxuid. The acoustic 
pressure and particle velocity can be four.d from a velocity potential 
solution using Eqs. (3.1.36) and (3.1.37). Consequently, the acoustic 
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pressure and particle velocity at the combustor exit can be determined 
from the velocity potential solution in the combustor. The acoustic 
pressure and particle velocity at the combustor exit can then be used 
to find the acoustic pressure and particle velocity at any other point 
in the ducted combustor system with four-pole equations in transmission 
matrix form derived using Eqs. (3.1.36)* (3.1.37). and (3.1.i»l). The 
combustor exit iiq)edance can be calculated from the duct exit impedance 
using the same foxu^pole transmission matrix equations. The transmis- 
sion matrix approach is discussed by Parrott (1973, Lampton (1978), 
Igarashi and Toyama (1958), Miwa and Igarashi (1959), and Igarashi and 
Arai (i960). 

Details of the mathematical models used in calculating the auto- 
spectra and cross-spectra are presented in the following sections. The 
next section discusses the source region. 

1. Source Region 

The acoustic state vector approach developed in the last section is 
used to study the spectral structure of pressure measurements made in 
the ducted combustion system. The linear non-homogeneous velocity po- 
tential wave equation assumed to apply in the source region is 

s G(u),x) ( 3.2.1) 


The boundary conditions at the combustor inlet and exit are, respec- 
tively , 


P c 
o o 


0 c 
o o 




x=0 


(3.2.2) 


and 
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CE 
o o 


O O 1 


x»L 


( 3 . 2 . 3 ) 


in the roiiovints iumJyuiu, ttu^ uumbuulioii nuisc uuurce is assumed to be 
acoustically compact in the x-direction. That is* the size of the 
source region is assumed small in comparison with the length of the 
waves generated. The volume source-distribution of monopole order con- 
fined to the duct cross-section x » is written as 

G(u,x) » C^JfXoj)6(x - *£)(l “ (3.2.1*) 

Similar source descriptions have been previously applied by Mani (1978) 
to a fan noise problem and by Morfey (1971 )» Ingard and Singhal (1975* 
and Swinbanks (1975) to duct noise source problems. 

Substituting Eqs. (3.1.38) and (3. 2. J») into Eq. (3.2.1), and the 
time Fourier transforms of Eqs. (3.1.36) and (3.1.37) into Eqs. (3.2.2) 
euid (3.2.3) and some algebredc manipulation yields a mathematical prob- 
lem that has the form of a general second-order Stuna-Liouville differ- 
ential equation with unmixed boundary conditions: 



(3.2.5) 


( 3 . 2 . 6 ) 


(3.2.7) 


where 
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c 



( 3 . 2 . 8 ) 


and 


e(-ik ) 

C = V (3.2.9) 

1 - >r 

s 

However, since the differential equation has coaq>lex coefficients the 
operator ^ is not self-adjoint or Hermitian. 

Since the differential equation is not self-adjoint, the develop- 
ment of a solution in terms of an orthogoned set of eigenfunctions is 
more conqjlex than obtaining a solution by constructing the Green's 
function solution (see Swinbanks (l9T5)» Morse and Feshback (1953), 
pp. 88U-886), and Kraft (1977) )• The Green's function solution is con- 
structed using the two initial value solutions of the homogeneous equa- 
tion 


= 0 ( 3 . 2 . 10 ) 

which satisfy Eqs. (3.2.6) and (3.2.7). This method is described hy 
Morse and Feshback (1956, pp. 523-530) and Friedman (1956) who apply it 
to self-adjoint operators. The solution to Eq. (3-2.10) with boundary 
condition fuven by Eq. (3.2.6) is 




* a 


(ik^O) 




(ik 0)' 
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X 


( 3 . 2 . 11 ) 


where 
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( 3 . 2 . 12 ) 


The solution to Bq. (3.2.10) with boundary condition given by Eq. (3.2.7) 
is 


I 

I 


I 


where 



(3.2.13) 


(3.2.1U) 


Consequently, the Green's function solution to Eqs. (3.2.5) to (3.2.7) 
is given by 




0 < X < X, < L 
i c 







0 < X < X < L 
I c 


(3.2.15) 


where W(x) is the Wronskian 


W(x^ = \l> — ^ ^ ^ 

' ^CI dx ''^CE dx 


(3.216) 
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Using the time Fourier trvisforiBs of Eqs. (3>1<36) and (3>1<37)» the 
pressure and particle velocity Fourier transforstt at the eoaibustor exit 
are 






(3.2.17) 


X*L 


and 

\ - ~ 5f (u.x.Xj^) ( 3 . 2 . 18 ) 

Equations (3.2.17) and (3.2. iP) define the acoustic state vector at the 
ccugibustor exit determined from the conibustor inlet and exit boundary 
conditions and the source spectrum. Relating the acoustic state vector 
at the combustor exit to the acoustic state vector at any other point 
outside the source region is done using an acoustic transmission matrix. 
The acoustic transmission matrix is based on the solution to the homo- 
geneous wave equation 'uid it will be discussed next. 

2. Duct Acoustic Transmission Matrix 

The duct acoustic transmission matrix relates the acoustic pres- 
sure and particle velocity at one point in the duct to the acoixstic 
pressure and particle velocity at another point providing neither point 
is in a source region and the duct has a constant area. It is calcur- 
lated as follows. 

Substituting Eq. (3. 1.^*1) into Eqs. (3.1.36) and (3.1.37) yields 
the pressure and particle velocity equations 

t (ik ii)*x (ik n)"jn 
^liae ° ♦ a^^be J (3.2.19) 

and 
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(ik) 

o 


ki*' 




.j2t» 


(ik^n) 




( 3 . 2 . 20 ) 


where 

■ 1 - M^R* ( 3 . 2 . 21 ) 

a^2 » 1 - (3.2.22) 

a^. » R* (3.2.23) 

and 

aj2 » n' (3.2.2l») 


The integration constants a and b in Eqs. (3.2.19) and (3.2.20) 
are determined by assisning that x ■ 0, p » p(0), and u ■ u(0) at the 
upstream bovmdary and x * L, p * p(L), and u = u(L) at the downstream 
boundary. As a consequence, the relation between the acoustic state 
vector at x * 0, y(0) and the acoustic state vector at L, Y(l) is 



where 


11 


(ik n)'*'L 

, ^ 11 *^ 22 ^ ~ '^ 12 *^ 21 ^ 


(ik q)'l 

o c 


® ll^'’2 " ^'’ l*12 


(3.2.26) 
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o,, ■ p c 

1 .; o 


(ikn)\ 

••a ae ^ ^♦aae ^ ^ 


® 11®22 ■ ® 21®12 


(3.2.27) 
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(ik fi)'*'L (ik a)~L 

oc 

21*22* ■ *22*21* 
‘’0*0^*11*22 “ *21*12 


(3.2.28) 


and 


22 


~* 21 * 12 * * 11 * 22 * 

*n*22 " *21*12 


(3.2.29) 


How that the transmission matrix between two points in a constant area 
duct has been derived, it is necessary to discuss the trcuismisslon ma- 
trix that applies across a discontinuity. 

3 . Discontinuity Transmission Matrix 

The four-pole trMsmission matrix used at the Junction of two pipes 
of different diameter is based on an acoustic energy conservation law. 
Various definitions for the acoustic energy in flowing fluids are given 
by Morfey (1971), ^fohri^g (1971), Candel (1975), Eversman (1979), 

Tester (1973), and Bergman (19^*6). While none of these definitions is 
applicable to a plane wave propagating through a soot particle cloud in 
a flowing field, the method for finding an energy conservation equation 
used by Morfey (1971 ) and by Bergman (19^<6) is applicable. The rela- 
tionship which has the form of an acoustic energy conservation law is 
derived using Eqs. (3.1.3**) and (3.1.35). A quantity corresponding to 
acoustic energy flux is defined by 
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N = 


(3.2.30) 


where the acoustic energy per unit mass is 


>1 ■ Vi * 


(3.2.31) 


and the mass flux fluctuation rate is 


•l “ "o”! * Vl 


(3.2.32) 


If Eq. (3.1.3^) is multiplied by u/c^ and Eq. (3.1.35) hy 

sum of the resulting second-order equations can be put in the form 


(-iu)) 
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s (-M = 0 

(ix Ip c I 
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(3.2.33) 


where E represents an acomtic energy density 
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(3.2. 3k) 


_ ^ 

pc o 
o o 



(3.2.35) 


Integrating Eq. (3.2.33) over a volume yields 


I 


(_im) (JV + 


f ^ n 

J s 


do = 0 


(3.2.36) 


since 


y*V-^dV= y* F-n 


do 


(3.2.27) 
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The area discontinuity can be assumed to take place sharply at the Junc- 
tion of two pipes with different diameters causing the flowing gas to 
undergo a sudden expansion or contraction. As a consequence of 
Eq. (3.2.36), across the pipe discontinuity the acoxistic energy is 
constant 

N(l)S(l) = N(2)S(?) (3.2.38) 

To derive a transmission matrix across the discontinuity Eq. (3.2.38) is 
rewrit t:.n as 

,/^d) * /^(2) (3.2.39) 

and 

S(l)w^(l) = S(2)w^(2) (3.2. 1*0) 

Algebraic majiipulation shows that the acoustic state vectors on either 
side of the discontinuity are related by 

Y(x 1 = [B ]Y(x^) (3.2. 1*1) 

1 Tn,n 2 

where the resulting transmission matrix is given by 


fB ] = 
‘ m,n 



(3. 2. 1*2) 


'( 2 ) 
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The next section discusses the last step in foraulating the model which 
is the specification of boundary conditions at the c^bustor inlet and 
the duct exit. 

System Boundary Conditions 

First, the boundary condition at the duct exit is discussed. Ex- 
perimental investi^tions of the exit impedauice of flow ducts reported 
by Mechel, Schilz, and Dietz (1965) and Cusmings (1978) show the duct 
exit pressure reflection factor may be greater them unity. Theoretical 
investigations reported by Lumsdaine (1977), Lumsdaine and Ragab (1977), 
and Mungur and Plumblee (1979) show that the pressxrre reflection factor 
is greater than unity became of the flow field at the duct exit. The 
following empirical expression for Ir| is given by Ingard and Singhal 
(1975) 


r| = 0.95 


(1 - M 

0 

* « 0 >. 


kr^ <0.5 


(3.2.U3) 


where |r| is measured at the upstream end of a duct. In the case of Jet 
exit flow, is replaced by (Lumsdaine (1977)). For the case 
studied herein |r| is about 1.03. 

The model calculationD shown are made using the following combus- 
tion duct exit pressure reflection factor 


F = 



- "q) 

- «o> 



(3.2. UU) 


the duct exit impedance is 
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M ♦ 



(3.2.l»5) 


Where 

V = 2kr^ j ^ ^ (3.2.U6) 

Jj^(v) is the Bessel function of the first order and first kindi S^(v) is 
the Struve function of the first kind and first order, and r^ is the 
duct radius. The duct exit impedance used was derived hy Lumsdaine and 
Ragab (1977) for a circular duct with flow having an open end fitted 
with an infinite acoustically rigid flange. Again, in the case of jet 
exit flow, is replaced by in Eqs. (3.2. UU) to (3.2.U6). 

The other boundary condition is specified at the combustor en- 
trMce. To model the impedance the combustor entrance is taken to be 
closed by a rigid circular plate at the entrance. This corresponds to a 
zero particle velocity at x = 0. Consequently, the pressure reflec- 
tion factor is unity at x = 0. 

Now that the elements that are needed to calculate pressure auto- 
spectra and cross-spectra have been described, the calculation proce- 
dure will be presented. 

C. Calculation Procedure 

The following procedure is used to calculate the pressure spec- 
trum. It is applied at each frequency as necessary to obtain the de- 
sired spectnsn. First, the exit pressure perturbation is arbitrarily 
assumed to be one Pascal. Then using the duct exit acoustic impedance 
based on the duct exit pressure reflection factor given by Eqs. (3.2.UU) 
to (3.?. ^6), the particle velocity i.s calculated. Next, using the duct 
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transsisBion matrix and the area discontinuity transmission matrix as 
necessary the acoustic state vector at the duct exit is used to find the 
acoustic state vector at the combustor exit. The resulting acoustic 
pressure and particle velocity are used to calculate the combustor exit 
impedance. 

The next phase uses this impedance, the combustor entrance impe- 
dance, and a white noise source spectrum to determine the velocity poten- 
tial solution in the source region. This velocity potentied solution 
is used to calculate the acoustic pressure and particle velocity at 
the combustor exit due to the specified source and boundary conditions. 

The last phase uses the duct transmission matrix and the area dis- 
continuity transmission matrix as necessary to calculate the acoustic 
state vector at any point in the duct from the acoustic state vector at 
the combustor exit. The pressure auto-spectra at a given point is cal- 
culated from 

Hii(u)) = ( p»(tij,x^)p^(u),x^)> (3.3.1) 

and the cross-spectrum between two points is calculated from 

= < p»(jj,x^)p^(w,Xj)> (3.3.2) 

The pressure level at a given frequency is calculated from 

PL^(oj) = 10 (o))l (3.3.3) 

In this chapter the acoustic model which will be used to calculate pres- 
sure auto and cross spectra has been described. The next chapter will 
describe the experimental investigation. 


Symbols Used in This Chapter 
A,B,C»D matrices 

velocity potential wave equation coefficients 
pressure and velocity wave equation coefficients 
boundary condition operator 
source spectrum coefficient 
sound propagation velocity, m/sec 
isentropic speed of sound, m/sec 
gas specific heat at constant pressure, J/kg-K 
diameter, m 

3 

acoustic energy density » J/m 
frequency, Hz 
solution in source region 
source term 

source term coefficient 


a 

^0 

c 

c 


d 

E 

f 

{(o,x,y) 

G(oj,x) 

H 


heat transfer coefficient, W/m -K 


J 

k 


L 

c 

M 

m 

% 

II 

n 


energy per unit mass, J/kg 

propagation wave nmber, ei ^ 

second order space differential equation operator 

combustor length 

Mach number, u /c 
o 

soot particle mass, kg 

heat transfer Nusselt number, Hd/ic 

acoustic energy flux 

number of particles per unit volume 


PL 


pressure level, dB 



X K X C ri- to ^ ^ ^ 
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P 

Q 

Q(-iw»ik tj) 
o 


'(u>)/c 


m,n 


m,n 


Y 

A 

e 

C 

n 

e 


pressure I N/m*" 

heat transferred to gas from particles by convection 

time and space Fourier transforti of q 

system state vector 

gas constant, J/Kg-K 

reflection factor 

radius, m 

entropy soxurce transfer function 
2 

area, m 

entropy of gas, J/kg-K 
ten?;erature , K 
velocity of bulk gas, m/sec 
Wronskian 

O 

mass flux, kg/m^^ sec 

Cartesian coordinate, m 

acoustic state vector 

acoustic impedance, mks Rayles 

duct transmission matrix element 

discontinuity transmission matrix element 

specific heat ratio of gas 

time and space Fourier transform of D/D6 

see Eq* (3*2. 8) 

see Eq. (3*2. 9) 

see Eq. (3*1 *18) 

time, sec 

gas thermal conductivity, W/m-K 
soot particle mass fraction, nm/p 
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p 

T 

S 

♦ [ 

♦ 

Q 


2 kr/(l - 

gas density, kg/n 

soot particle thermal relaxation time, 

Fourier transform operator 
reflection factor phase euigle, deg 
velocity potential function 
velocity potential wave number factor 
angular frequency, radians/sec 


Superscripts and subscripts: 


< > 
C) 

n 
( )■ 

( )" 

{ )° 

^ ^CE 
( ^CI 
( >d 

( ). 


( ) 

r 

{ ) 

c 

( 

( ) 


n,n 


time average 

vector quantity 

instantaneous quantity 

upstream propagating 

downstream propagating 

wave propagating with flow 

combustor exit 

combustor inlet 

exit of duct system 

identifies an axial duct location 

identifies source location 

elements of a matrix 

property of or due to soot particle 

reference state quantity 

perturbed quantity 

property far from soot particle 


W 



I 


r 
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IV. 13CPERIME3ITAL IHVESTIOATION 
A. Apparatus, Test Conditions, and Instriraentation 

The test facility is shown schematicalxy in Figure 5. The conn 
hustor consists of a J-l*7 burner can placed concent ricedly in a 0.30 m 
diameter by 0.77 m long flow duct. The combustor section is followed 
by a 0.38 m diameter by 0.76 m long spool piece. This section is fol- 
lowed by a 0.30 m diameter by 6.1 m long flow duct. A photograph of 
the test facility and test site is shown in Figure 6. The nozzle 
shown mounted to the long duct in Figiu*e 6 was removed for the tests. 

The measurements discussed were made at an exit temperature of 
approximately 9^’0 K and at air mass flow rates of 0.5» 1*13, and 1.68 
kg/sec. The corresponding air velocities at the exit of the long duct 
were l8.5> **1.6, and 6l.3 m/sec and the corresponding fuel flow rates 
were 0.009* 0.018, and 0.027 kg/sec. The fuel-air ratio was about 
0.02 for each test condition. Also, at each of these test conditions, 
the fuel used was a Jet aircraft fuel called Jet A. 

Sinraltsuieous internal fluctuating pressure measurements were made 
at the three locations shown in Figure 5. One measurement was made in 
the spool piece. Two measurements were made in the long duct, one Just 
downstream of the spool piece and one near the nozzle exit plane. 

The transducers used were conventi<->nal 5/8 cm diameter (nominal) 
condenser microphones with pressure response cartridges. To preclude 
exposing the microphones to the high temperature combustion gases in 
the duct system, the microphones were mounted outside the duct. The 
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fluctuating static pressure at a given location in the duct system vas 
connected to a microphone by means of a "semi-infinite” acoustic wave- 
guide tube which was a very long piece of narrow tubing. This design 
eliminated longitudinal resonance and produces am anechoic connection 
between the duct system and the transducers. 

A drawing of a typical probe is shown in Figure 7. The micro- 
phone was flush mounted in the acoustic waveguide through a supporting 
block euid housed in a pressure chamber. Attached to the block were a 
5/8 cm diamieter sensing tube on one end and a coil of tubing of the 
same diameter, 30 meters long, on the other. The diameter of the coil 
was approximately 30 cm. 

The sensing tube of each probe was flush mounted as a static 
pressure tap at each measuring location. A regxilated nitrogen purge 
flow was maintained in each sending line to protect the microphone 
from hot combustion gases, fiincc the signals frcm pairs of probes were 
to be used to calculate cross-spectra, the regulated purge flow system 
was separate for each probe. This prevented any common valve noise 
from the regulator contributing to the cross-spectra. Static pres- 
sure was balanced across the microphone by means of a small vent hole 
connecting the pressure chamber and sensing tube. A schematic of a 
typical duct probe installation ir, shown in Figure 8. A photograph of 
duct probes installed in the spool piece is shown in Figure 9* 

The internal probes have previously been used for engine measure- 
ments (Karchmer and Reshotko (1976), Peshotko, Karcluner, Penko, and 
McArdle (1977), and Karchmer, Reshotko, and Montegani (l9T7) and mea- 
surements in a combustion component tent facility (Reshotko and 
Karchmer (1977). Probe design, frequency response , and operating char- 
acteristics ore discussed by Karchmer (1978). 










78 


B. Data Acquisition and Processing 

As part of the facility start up procedure* the microphones vere 
calibrated with a pistonphone which provided a standard sound pres- 
sure level of 12I4 dB (referenced to 20 yPa) at 250 Hz. Microphone 
signals were sent over shielded cable to a remote control room for am- 
plification, signal conditioning, and recording on magnetic tape. 
The tape recorder used was equipped with IRIG standard, intermediate 
band, group II electronics, which provides a flat frequency response 
to 20 kHz on both record and reproduce at a tape speed of 60 in. /sec. 
The dynamic range available was approximately ^8 dB. All signals were 
recorded on 120 second record lengths for later processing. Record 
lengths of 120 seconds were used to obtain high statistical confidence 
in the spectral data. Spectrxam function confidence limits are de- 
scribed in Appendix B. All measurement errors are determined by the 
microphone response which is generally accurate to better than plus or 
minus 1 dB. / 

Operating conditions were monitored suid recorded on the NASA 
Levis CADDE system (Central Automatic Digital Data Encoder). On line 
data processing and analysis was limited to 1/3 octave-band spectral 
analysis of the various pressure signals. 

All auto-spectra and cross-spectra produced from pressure signals 
were obtained by off-line processing of the taped signals on a twc- 
channel, FFT (Fast Fourier Transform) digital signal processor, witi" 
built-in aneLlog to digital converter and 120 dB/octave anti-aliasing 
filters. Tiie 99 percent confidence interval for the spectra computed 
from the taped signals is shown in Appendix B to be better than plus 


or rainus 1.5 dB. 
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The auto-spectra and cross-spectra considered herein are defined 
as follows. The finite Fourier transform of time signals x(B) and 
y(9) of the record length 0 is 



( 1 *. 2 . 1 ) 


(U.2.2) 


The two one-sided auto-spectral density functions (also called power- 
spectral density functions) and the single one-sided cross-spectral 
density function obtainable from measurement of signals x(0) and 
y(6) are defined by the following expressions 


G,^(f) = 

lim 2E 

[X»(f)Xjj(f)] 

(i*.2.3) 

|o 

8 

t 

0 

G^(f) . 

lim 2F. 


(h.2.h) 

0 -► w 0 


lim 2E 

• [X»(f)Y^(f)] 


0 -► 00 0 



= C (f) + 
xy 

iQ 

xy 

(i*.2.5) 


where E[ j denotes nn expected value operation (Bendat and Piersol 
(1971)). The auto-spectra arc shown in the following figures in terms 
of pressure levels referenced to PO uPa defined as 

= 10 logjQlG^(f)) (U.2.6) 

PLy = 10 log^Q[Gy(f)J 


(U.2.7) 
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The cross-spectra are shown in the following figures in terms of pres- 
sure levels referenced to 20 yPa and phase mgles defined as 

PI,^(f) - 10 108, 



(f) * tan 


-1 








C* Ffeasured Pressure Auto-Spectra and Cross-Si>ectra 

Constant bandwidth pressure auto-spectra and cross-spectra mea- 
sured in the spool piece, near the entremce of the long duct, and necir 
the exit of the long duct are shown in Figures 10 to 13. The spectra 
at each location have a different structure. However, at each loca- 
tion, the structure of the measured spectra is similar for each test 
condition. For example, the location of resonance peaks and dips in 
Figures 10 to 12 is nearly the seune for each operating condition. At 
each location and for all test conditions the peaks and dips tend to 
be sharper at the low frequencies and broader at the high frequencies. 

Constant bandwidth pressure cross-spectra measured across the 
area contraction at three test conditions are shown in Figures 13(a), 
(b), and (c), respectiveJy . Both magnitude and phase arc shown. The 
location of the resonance peaks and dips in th'^ magnitude plot of the 
cross-spectra is nearly the same at each test condition. Also, the 
peaks and dips again tend to be sharpter at the low frequencies and 
braoder at the higher frequencies. 

The phase plot can be characterized as neither Just a negative 
sloped straight line due to a time delay nor Just a curve that varies 
abxniptly between 0 and l80 degrees in a regular manner that would be 
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Figure 11 - Pressure auto-spectra measured near entrance o( long 
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V. 
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Figure 12. - Pressure auto-spectra measured near exit of long duct. 
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due to the presence of waves havinc approximately equal ma^itudes but 
opposite sign travelling upstream and downstream. Instead the phase 
plots show both tendencies with the time delay tendency dominating at 
the higher frequencies. 

A large peak appears above 200 Hz in the auto-spectra and cross- 
spectra measured at air mass flows of 1.13 and 1.68 kg/sec. These 
peaks appear at such high frequencies that they cannot be related to 
specific duct resonance harmonics. These peaks may be due to a feed- 
back resonance between the combustor and the duct. They are not ana- 
lyzed herein. 

In the next chapter calculated auto-spec '^ra and cross-spectra will 
be compared with the measured auto-spectra and cross-spectra discussed 
in this chapter. 

Symbols Used in This Chapter 


xy 
E[ ] 


xy 

i 

PL 

*^xy 

Xo(f) 

x(e) 

y(e) 


real part of cross-spectra 
expected value operator 
frequency 

auto-spectra of x signal 
auto-spectra of y signal 
cross-spectra between x and y signals 

pressure level 

imaginary part of cross-spectra 

finite Fourier transform of x time signal 

X tim<‘ signal 

finite Fourier transform of y time signal 
y time signal 
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0 record length* aec 

0 time, nee 


xy 


croHD-niieetrum phano amUr 


V. C30MPARIS0M WITH EXPERIMEHTAL DATA 


The single time constant model discussed in Chapter III is used 
now to calcxilate auto-spectra and cross-spectra at locatiuris which cor- 
respond to the locations of the measured auto-spectra and cross-spectra 
presented in Chapter IV. The parameters used to calculate the mass 
fraction and soot particle relaxation time used are given in Table III. 
The value of the soot particle relaxation time is considered to be 

1.5*<10 ^ sec for all cases. The time constant, x /tc , used in each ex- 

s s 

ample then is due to the presence of different soot particle mass frac- 
tions as shown in Table III. 

For a wave propagating downstream suid for the time constants given 
in Table III, the sound propagation speed and attenuation calculated by 
substituting Eq. (3.1.25) into Eq. (2.1.t»l<) and Eq. (2.1.U6) are iuown 
in Figure 1^4. For a time constant of 1.0 sec, the sound propagation 
speed is 6l0 m/sec which is the isentropic sound propagation speed. 

The sound propagation speed is above 600 m/sec above 50 Hz for a time 
constant of 0.010 sec. However, for a time constant of 0.001 sec, the 
sound propagation speed varies from the isothermal speed of sound prop- 
agation, 510 m/sec, to the isentropic speed of sound at high frequen- 
cies. For a time constant of 0.0001 sec the sound propagation speed is 
the isothermal speed of sound below 3^0 Hz. 

For time constants of 1.0, 0.010, and 0.0001, the attenuation is 
less than 0.5 dB/m at most frequencies below 3^0 Hz. However, for a 
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TABLE III. - Parasieters used to calculate dispersion and attenuation 
due to soot particle-gas heat transfer 


t, K 


K, W/o-k 

. . . . 5.38x10“^ 

Pq. kg/M^ 

0.378 

Cp, J/kg-K . . . . , 


y ........ . 

i.i» 

0 . kg/M^ 

s 

1 ■''n.n 

^s» 

. . . . 0.3^2x10“^ 

c , lii/sec 

o 


“s* 

. . . . 3.15x10”^^ 





Example 



1 

Isothermal 

2 

Mixed 

3 

Adiabatic 

k 

Adiabatic 

n, number /m 

1.8x10^^ 

1.8x10^^ 

1.8x10^^ 

1.8x10^^ 

nnig, gm/m 

5.67 

0.567 

0.0567 

5.67x10"** 

K , Nm /p j 

s s 0 

1.5x10”^ 

-U 

1.5><10"^ 

ixio"^ 

1.5*10 

-2 

1.5x10"^ 

T /< , sec 
s s 

1x10 

i 

IxlO 

1.0 












FREQUENCY. HZ 

Figure 14. - Cilculated attenuation and dispersion. 
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tine conft«nt of O.OOl sec» the attenuation is greater than 1 dB/a 
over the freqvrency range of 150 to 3U0 Hz. 

A. Auto-Spectra 

Calculated auto-spectra are con^ared with auto-spectra Beas\u*ed 

with a mass flow rate of 0.5 kg/sec in Figures 15 to l8. Auto-spectra 

calculated with t„/k^ ■ 1.0 sec and 0.001 sec for the location in the 
s s 

spool piece upstream of the area contraction are compared with auto- 
spectra measured at that location in Figures 15(a) and (b). For 

■ 1.0 see, the sound propagation speed is isentropic and the at- 
tenuation is negligible. Consequently, the dashed curve shown in Fig- 
ure 15(a) presents the result which would normally be predicted. The 
calculated dips and peaks shown in Figure 15(a) are in different loca- 
tions than the ones measured. Also, the dips an. peaks are harp at 
all frequencies while the measured dips and peaks are sharp only at the 
low frequencies. 

However, for * 0.001 r.ec, the calculated dips and peaks are 

s. s 

in better agreement with the measured dips and peaks. In addition, the 
dips and peaks are sharp at low frequencies and broad at high frequen- 
cies in excellent agreement with the measured spectre. The improve- 
ment in the position of the calculated dips and peaks is due to the 
sound propagation speed calculated with the time constant of 0.001 sec. 
Thir, time const.ju)! pmiuce:^. the nonlinear beh.avior of the sound propa- 
ftation npeed wlii ‘h the mito-rpoo t ra ourve ohovn in Firure 

in a n<''n--im-ar fashi(*^n. Furthermore, the ‘hange in time con- 
^^tanls fr ra to 0.001 ehangeu lh*‘ ral i*tructure from hav- 
ing iharp -u i dip.: at nil i*r<^queri'*l having sharp peaks and 

dips at fr” i and br'*ad arid ilpr 'U high f requeneien . 
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Figure 16. - Comparison of measured and calculated auto-spectra 
near long duct entrance (Ug • 18.3 mfsec). 
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Figure 17. • Cwnparlson ol mMsurcd and calculaled pressure spectra near 
entrance to long duct (Uq • 18.5 m/secl. 
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IfeMxired auto-spectr« at the location near the long duct entrance 
dovnatrean of the area contraction are coatpared vith auto-apectra cal- 
culated using tine constants of 1.0 and O.CKA sec in Figures l6(a) 
and (b), respectively. Again, calculati<ms aaide using the 0.001 aec 
are in better agre«Beot with the Masured data. 

1%at \ising the isotheraal speed of sovmd propagation <«ther than 
the isentropic speed of sound propagation in model calculations does 
not produce auto-spectra %^ich agree with masured auto-si>ectra is 
shown in Figures 17 and l8. In these figures measured and calculated 
spectra near the entrance and near the exit of the long duet are c<os> 
pared. Auto-spectra calculated using a time constant of 0.0001 sec, 
which produces little attenuation and an isothermal sound propagation 
speed below 3^0 Hz, are shown in Figures 17(c) and l8(c). The calcu- 
lated curves are cos^ressed versions of the curves shown in Figures 
17(a) and l8(a) which are calculated with a time constant of 0.01 sec 
which produces little attenuation and an isentropic sound propagation 
speed above 80 Hz. Neither the curves calculated with a time constant 
of 0.01 sec or 0.0001 sec are in as good agreement with the measured 
data as the curves calculated using a time constant of 0.001 sec shown 
in Figures 17(b) and 18(b). The changes with mass fraction appear 
graphically in the form of stretching or as a compression of the basic 
structure. These changes were easy to recognize but proved difficult 
to simply explain. 

In summary, the single time constant model developed in Chapter III 
provides a simple, elegant theory which explains the spectral structure 
of the measured auto-spectra. In addition, the next section shows it 
tQso explains the structure of the measured cross-spectra. 
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B. Crosa-Spectx^ 

Ibe cz^BSo'tpttctra calculated using tiae constants of 1*0 and 
0.001 sec are compared with t)w data Measured for « l8 a/see In 
Figure 19. ^n>e aodel calculations aade using a tias constant of 
1.0 sec shovn in Figure 19(a) produces a cross-spectra vbicb does not 
look like the aeasured cross^spectra. However, the cross^spectra cal- 
culated using a tine c<mstant of 0.001 sec shovn in Figure 19(h) does 
res^hle t)M aeasured cross- spectra. The peaks and dips in the cross- 
spectra aa^iitude calculated using a time constant of 0.001 see are 
close to the locations of the peaks and dips in the aeasxn'ed cross- 
spectra. In addition, the structure of the measured and calculated 
magnitude curves is similar. Furthermore, the phase angle calculated 
using a time constant of 0.001 sec and shovn in Figure 19(h) is in ex- 
cellent agreement with the structure of the measured phase angle. 

C. Ifean Flow 

Calculations made using a time constant of 0.001 sec with three 
different mean flow velocities are shown in Figxire 20. The calculated 
pressure level at the exit of the long duct changes slightly with the 
meari flow speed at the low frequencies. Vfhile the change is greater 
at the higher frequencies it occurs predominantly in the spectnm 
level and not in the frequency location of the peedcs. Consequently, 
these changes could not be observed in the nweisured spectra. However, 
the measiu*ed spectra shown in Figure 11 for the three flow velocities 
have similar trends at low frequencies. Accordingly, the spectra mea- 
sured at the higher duct exit flow velocities are also in fair agree- 
ment with the model calculations made for the time constant of 


0.001 sec. 


mmmm 
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EXPCRIMENT 



Figure 19. - Comparison ol measured and calculated cross- 
spectrum across area contraction (Uq • 18. 5 m;sec). 


PRESSURE LEVEL dB 



FREQUENCY. Hz 

Figure 20. - Effect of flow velocity on comp^ed pressure 
spectra: TjIkj -0001 sec, Kj • L5xl0 ^ 
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D. Discussion 

No attempt vas made in this investigation to detemine soot part- 
icle size or nass concentration. Data In the literature indicate that 

typical soot particle diameters range from 0.05 to 1 vm (Linden and 

3 

Heywood (1971)). Smoke concentration ran^s fran 0.53 to 5*8 gm/m 
in the primary combustion zone and from 0.00018 to 0.015 gm/m'^ at a 
combustor exhaust station were found by Norgren (1971 ). Besides being 
a function of position in the combustor, smoke concentration was found 
to be a function of combustor model, fUel-air ratio, and operating 
pressure by Norgren (l97l). 

1. Spectra Structure 

The spectral structure for a given operating temperature and a 
given geometry is determined by the value of the ratio of the soot par- 
ticle thermal relaxation time to the mass fraction. Using Eqs. (2.3.9) 
and (2.3.12) and the parameters in Table III, the soot particle concen- 
tration can be calculated from the soot particle radius for a given 
value of this time constant by 

nm * l».81*xio^(r^V(x /< ), gm/m^ (5.1) 

s \ s// s s 

This function is plotted in Figure 21 for three values of the time con- 
stemt. Also shown in Figure 21 is the smoke number for a given soot 
particle concentration as given by Norgren (1971 ). The circled point 

in Figure 21 at r = 0. 3**2 um and (nm ) = 0.567 gm/m corresponds to 

s s 

the value in Table III for x^/tc = 0.001 sec which produced the spec- 

•J s 

tral structure that most resembled the measxired structure. This par- 
ticular point may be unrealistic since it corresponds to a large 

smoke number. However, any point on the x /ic a 0.001 sec line woxild 

s s 
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Figure 21. - Parameter map for Tj/Kj. 
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alto prodttoe the taae tpectrel struotxare. Contequmtly, theee retultt* 
to this extent* ere independent of the pareiMter in fthle III. Unfor- 
tunately* this also Beans that the aethod cannot Independently he used 
to eiitiaate the soot particle radixis or mass concentration. 

Figure 21 also indicates how the sound propagatim speed raries 
with soot particle radius and mass concentration at ccnhustion noise 
freq\Mncies. The sound propagation speed can he detendned frca Eqs. 
(2.1. and (3.1.25) to (3.1*28). The sound propagation speed is iso- 
theraal at frequencies for which (mt./k ) is less than unity. Conse- 
quently* sound propagates isothemally at frequencies less than a cor- 
ner or break frequency* f* given by 

8 8 

Above the comer or break frequency the sound propagation speed is 
changing from isothermal to isentropic. Along curves A, B* and C 
the break frequencies are, respectively* 15.91* 159.1* and 1591 H*. 

For a time constant of 0.01 sec, corresponding to curve A in Fig- 
ure 21, the speed of sound propagation exhibited in Figure lU is less 
than 5^0 m/sec at frequencies below 15*91 Hz. Also, for a time con- 
stant of 0.001 sec, corresponding to cxunre B in Figure 21, the speed 
of sound propagation exhibited in Figure iJ* is less than 5^0 m/sec at 
frequencies below 159.1 Hz. Consequently, at comb\>stion noise frequen- 
cies, in the region above curve A in Figure 21 the sound propagation 
speed calculated is isentropic. In the region below curve C the cal- 
culated speed of sound is isothermal. However, along curve B a 
mixed speed of sound propagation occurs at combustion noise frequen- 


cies. 
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The aaximuB enount of attenuation is difficult to detemine froa 
Gq. (2. 1.U6). However, it is possible to determine the frequency 
which maximizes the phs«e angle of (ikO^). The attenuation at this 
frequency is about 5^ percent of the maximal attenuation. Substituting 
Eq. (3.1.23) into Eq. (2. 1.U6) and maximizing the phase angle shows 
the phase angle aaximuB attenuation at any operating tesq>erature , for 
M ■ 0.0, occurs at 


“ 2.{T /k ) 

8 8 


(5.3) 


the phase angle maximum attenuation at this frequency is 

o s I4.96 <iB/m (5.**) 

oiajc c 

o 

Thus, along curves A, B, and C in Figure 21, the phase angle maxl- 
mun attenuation occurs, respectively, at 18,83, 188.3, and 1883 Hz and 
is, respectively, 0.15, 1*53, and 15«3 dB/m. The attenuation is, re- 
spectively, O.28U6, 2.8U6, and 28. U6 dB/m at f » 1.0*10^^. 

The measured spectra shown in Figures 10, 11, and 12 may have 
similar structure at each operating condition because the value of the 
ratio of the heat transfer time constant to mass raction did not 
change much with operating condition and was approximately 0.001 sec. 

On the other hand, if the broad peak above 200 Hz, which occurs 
in the measured spectra taken with long duct flow velocities of UI.6 
and 61.3 m/sec, is due to feedback between the duct acoustics and the 
combustion process, then its occurrence at these operating conditions 
may be due to a decrease in acoustic damping of the longitudinal waves 
due to a change in soot mass fraction or radius at these operating 
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conditions. This possibility suggests th«t soot pa^tleles era be used 
to control certain types of coabustor instability ty attenuating pres* 
sure oscillations at critical fre<iuencles. 

2. Model Assumptions 

In Chapter III the soot particles are assumed to hare a unifon 
tesq;>erature and no mass loss because these simple assumptions produced 
theoretical resxilts which agreed with measurements, because they 
greatly simplified the analysis and because more conplicated assvnp* 
tions did not produce better models. The actual effect of a pressure 
perturbation on the teiq>erature and mass equilibrium between a liquid 
or solid particle and a geis is a complex problem depending on at least 
the following influences: 

1. Convective heat transfer to or from the particle. 

2. Heat transfer within the particle. 

3. Mass transfer to or from the particle. 

Surface chemical reactions on a soot particle. 

5. Liquid-vapor phase transitions of a droplet. 

6. Radiative heat transfer. 

Radiative heat transfer la not an important factor in this case because 
the combustion duct operating temperatures are low. 

Many models can be constructed which describe this problem for a 
soot particle. As an example, it is possible to analyze this problem 
by assximing instantaneous heat transfer within the particle and by 
using a pressure and temperature dependent surface chemical reaction 
which effect heat and mass transfer as was done in Chapter II. Tlii; 
model is complex and dependa on many parameters. 


10l( 

In Chapter III the particle* are asstawd to have a unlfora tem- 
perature. This 1* aasuaed to be due to some unspecified pressure and 
tesq;>erature dependent surface chemical reaction. The aodel developed 
is dependent on only a single parameter, the ratio of to 

It is possible to analyse this problem "by ignoring surface ctuni- 
cal reaction effects on heat and mass t..ansfer and by assuming ln« 
stutaneous heat transfer vithin the particle. Ibe resulting thermal 
relajmtion time for a l-pm particle is about 0.3 vcec which shows that 
this type of particle quickly follows any teiq>erature change and does 
not have a constant time independent tesq>erature. However, experl> 
ments show that smoke concentrations of soot particles eure reduced as 
soot travels from a measuring station at the primary cembustion zone 
to one at the combustor exhaust (Norgren (1971)). Consequently, sur- 
face chemical reactions cannot be ignored. Further study is needed to 
discover what assxsnptions are best for analyzing this problem. 

Symbols Used in This Chapter 

f. corner or break frequency, Hz 

b 

f frequency of phase angle maximum, Hz 

a phase angle maximum attenuation, dB/m 

m&x 
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VI. SlDMOff AHD CGHCUUmiG REIMISS 

A. Suoiuzy 

In Chapter I, the motivation for the ccwbined experimental and 
analytical p3rograB described in this dissertation is establisted. Spe- 
cifically, the importance of combustion noise in future aircraft design 
is indicated. Current approaches to combustion noise research are dis- 
cussed. Next, the need to study the propagation of combustion noise is 
established. Then the topics investi^ted herein are specified as be- 
ing concerned with cuisvering the following two questions: 

(1) What is the wave equation in an oi>erating ducted combustion 
system? 

(2) What formulation of the solution to the wave equation provides 
an acceptable mathematical model of an operating ducted coobustion sys- 
tem? 

Finsdly, the available work of iu>.:;rest in these eiresis is discussed. 

In Chapter II, the plane wave equation as modified for the re- 
si>onse of a fuel droplet cloud or a soot particle cloud is discussed. 

It is shown that at low frequencies sound can propagate at speeds near 
the isothermal speed of sound. 

In Chapter III, an acoustic model of the combustion duct system 
which depends on a single time constant is presented. This chapter 
describes the use of a two-port transmission line using transfer ma- 
trices to represent sound propagation in the homogeneous differential 
equation solution region of £ui operating ducted ccmibustion system. The 


■ 1 

I 

- t 

f 't 

i -I 

1 4 




i 


I 


i 

i 


105 


106 


Mae of an approach based on construction of the one-diisenaional Green's 
function In the non-hotaogeneous region Is discussed. 

Ctepter IV first describes the experlnental hardware* instruaeii- 
tation, and fre<iuency doswdn representation. Next* measured auto- 
spectra and cross-spectra are presented. Briefly* the ducted combus- 
tion system consisted of a J-W can combustor* followed by a short 
spool piece, fluid then by a long duct. Fluctuating pressure measure- 
ments within the operating ducted combustion system were made with 
speciedly constructed "semi-infinite*’ waveguide probes, with conven- 
tional microphones being used as the pressure treuisducers. The data 
obtained from these probes were off-line processed on a fast Fourier 
transform digltflO. signed processor. The processor wais used to obtfldn 
pressure level sjiectra as well as cross-spectra between internal pres- 
sure meamurements. 

In Chapter V the measured auto-spectra and cross-spectra eure com- 
pared with the calculated auto-spectra and cross-spectra and the re- 
sults are discussed. Auto-spectra and cross-spectra calculated using 
a time constant of 0.001 sec in the model developed in Chapter III are 
in excellent agreement with the measured auto-spectra and cross- 
spectra. 

B. Recommendations for Future Work 

The experimental data obtained in the present study support the 
idea that combustor emissions in the duct act as the primary mechan- 
ism producing the attenuation and dispersion of combustion noise in 
this operating liquid fuel cembustion system. However, combined 
acoustic and combustion measurements need to be made if the theory is 
to be compared with the experimental data. Thus, this Is an area 
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where continuing research efforts shovdd be directed. Unfortunately, 
answers to the irany reaaining questioiis will probably be found only 
through aultidisciplinary research. 

C. Concluding Remaurks 

The usefulness of a single paraseter linear systos aodel to de- 
scribe the dynandc behavior of an operating ducted liquid fuel c(»bus- 
tion system has been established. Vfhen allowance is made for disper- 
sion and attenuation effects and the boundary conditions are properly 
specified good agreei^nt between measured and autalytical spectra and 
cross-spectra is achieved. However, it is not clear if these results 
can be related or applied to actual engines. For example, the design 
of liners or suppressors for combustion noise reduction obviously re- 
quires knowledge of the appropriate wave equation arid the correct prop- 
agation speed of sound. Clearly, if dispersion and attenuation due to 
soot particles or liquid fuel droplets occurs in actual engines at low 
frequencies it must be tedcen into account. However, the propagation 
speed of sound at low frequencies in actual engines is unknown. 

Finally, as stated in the INTROIXJCTION , the objective was to use a 
dynamic system state-space approach to study the operating ducted com- 
bustion system. The measurement techniques described herein, as well 
as these mathematical methods have been available for many years. It. 
is the bringing together of both In a unified research program in a new 
area with practical objectives that lead to this work making a major 
contribution to the study of combustion noise. 
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RSSP01I8E OF THE FUEL DROPLET CLOUD 

The purpose of this appendix is to discuss the derivation of the 
transfer functions «4T(u)/p and JM/c for the case of sound 

w ^ 

propagation in a liquid fuel droplet cloud following the approach used 
hy Marble and Candel (1975) for sound propagation in a water droplet 
cloud. 

The following assumptions are made: 

1. The ft»l drop is spherically symnetric. 

2 . Convection corrections can be obtained with experimental cor- 
relations . 

3. Fluid properties are assumed constant with the property values 
determined at an appropriate reference condition. The liquid fuel is 
assumed to be a single pure component. 

U. The relationship between the gas phase fuel concentration and 
liquid fuel droplet temperature is given by the vapor pressure corre- 
lation for the pure liquid, i.e., the Clausius-Clapeyron equation. 

5. The ambient gases have negligible solubility in the liquid 
phase and only the fuel vapor is diffusing from the siurface. 

6. The radial motion of the liquid surface is assumed to be 
small. 

7. Heat transfer by radiation is neglected. 
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8. Mass diffusion is represented by an effective binary diffusion 
law corrected for convection* and the Dufour effect is neglected in the 
heat flux (Bird* Stewart* and Lightfoot (I960)* and Velty* Wicks* and 
Wilson (1969) )• The diffusion coefficient is calculated fro* 



and the diffusion Ifusselt number is given by 

HUjj - 2 + 6 Re^'^%c^'^^ (A-2) 

9* Reaction effects such as fuel oxidation are negligible. 

10. Heat transfer is calcinated using the Chapman>Ensog fonsula for 
the thermal conductivity of a monatomic gas at low density at tempera- 
ture t corrected for convection (Bird, Stewart, and Li(ditfoot (I 96 O)* 
and Welty, Wicks, and Wilson (1969))- The thermal conductivity is ced- 
culated frcai 

ic = 8.322Uxio"^ ^^*3) 

and the heat transfer Nusselt number is given by 

Nujj » 2 + 6 Re^'^^Pr^'^^ (A-U) 

Surface tension and siirface charge corrections are omitted. 

11. The fluid in the drops is rapidly mixing creating droplets of 
uniform temperature. 

12. The mean flow velocity is neglected in the vapor and bulk gas 


continuity equations. 
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A perturtetlon in pressure causes a perturbation in the hulk gas 
tenperature* droplet t^perature, the mass source » 

♦, I and the partial pressure of the fuel far fra» the drop, p 
The following system of equations is derived to relate these perturba- 
tions. This common factor is dropped henceforth. 

The bulk gas energy equation will be considered first. Since the 
heat conducted from n liquid fuel droplets of radius r^ into a unit 
voluae of gas having thermal conductivity k by a convective heat 
transfer process with Nusselt number Nu^^ is equeil to the product of 
teiq>eratui*e and the rate of entropy rise, the bulk energy equation 
is 



where the fuel droplet mass loading is given by nB^, the fuel droplet 
mass fraction, k^, is 

<L = (A-7) 

and the fuel droplet thermal relaxation time, is 


Tj^ = mj^Cpy' [l»irrj^ic(Nu^/2)] (A-8) 

Perturbations in bulk geis temperattire and pressure are related to 
perturbations of fuel droplet temperature next. A small perturbation 
expansion of the perfect gas state equation (Eq. (2.1.5)) yields 


Ill 


Pl/Po • Pi/Pq - (A-9) 

Equation (A-9) is then substituted into the sero mean flov, ■■>*11 per- 
turbation form of the ideal gas entropy equation (Eq. (2.1.8)) to 
eliminate the density perturbation. Finally, the tin» derivative is 
taken and the bulk gas energy equation (Eq. (A-6)) is used with the 
resiilt that 


A. 


de 

1^0 



(A-10) 


or in terms of the liquid mass fraction and thermal time delay 



- (y - 1) 


d Al\ 

'“*VV ' ‘o 


(A-11) 


The discussion so far has neglected the mass source term, The 
first step in intiroducing this term is to relate perturbations in vapor 
pressure, fuel droplet temperature, and mass flux to pressure pertur- 
bations. The time rate of change of the droplet mass is governed by 
the rate at which the fuel vapor can diffuse into the background gas. 
Thus, 



(A-12) 


The mass of vapor entering a unit volume occupied by n droplets is 
in terms of the partial pressure 


t = nl»itr,p D„ 

L o Vg 


%Pv-L - Py, 


> Py-L * 


‘'y-L 


(A-13) 
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For a small vapor pressure perturbation 


P - Po ♦ Pj 


a small mass addition occurs 


♦ ■ *0 * *1 


vhere 


*1 “ "^’'*'L%®vg\-r 


"d\ f Vl.1 ~ 


The Clausius- Cl apeyron equation is used to relate the fuel vapor 
pressure perturbation to the fuel droplet temperature perturbation as 
follovs ; 


S-L * Vl,1 


^-L.l 




since 


l/p » 1/p- 
CL 1 


^v.o * ‘’v 


Substituting Eqs. (A-17) and (A-l8) into Eq. (A-l6) yields 


♦i * V V— I ~an— 


i 
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lfov« vlth th« following themodynoBic relations: 


Cp - cy • «/(MM)g 


Y • Cp/c^ 




g®P 


Eq. (A-31) can be written as 


o \ 0/ 


J<wr-p 


VV2 ' 


Xo 

P. 


K-l-V V \(’*>y Vi h 
[vp/ - y^s ‘0 ■ ' 


or as 


»o ’d T’-l" ‘o ^,cj 


where the fuel droplet ness fraction, iCj^, is again given by 


* nm/p^ 
L o 


the molecular weight ratio of vapor-to-b\ilK gas, 6, is given by 


(MW), 


TmwT 


the latent heat parameter, n, is given by 


t c 
o p 


and the diffusion relaxation time, is given by 


(A-22) 

(A-23) 

(A-2M 



(A-25) 


(A-26) 


(A-27) 


(A-28) 


(A-29) 
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nie fuel droplet energy eqiwticm involves the mass flux perturba- 
tion as follows. The tiae rate of change of a fuel droplets energy is 
equal to the instantaneous difference between the convective heat input 
and the heat required for surface vaporization. Thus , 



or 


■ 


(A-32) 


where 


T 


T,L 



(A- 33) 


The fuel vapor mass source perturbation is related to the fuel 
partial pressure perturbation and the total pressure perturbation as 
follows. From the vapor contin\iity equation 


36 ^V,l 


+ 


3x 


(A-3U) 
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•Dd from the bulk gM contiouity equation 


d 

30 ^1 


♦ 


O 3x 


♦l 


(A- 35) 


Using Eqt. (A-3^) and (A-35)» the velocity gradient can be eliadnated 
as follows: 


3u. 


*1 1 

^ 'v.o ' "v.o ” 


P- ” P_ 3® 
o o 


(A-36) 


Thus, 


r 1 jj| , / 1 y*^v.i / Affi 

[v.o’^j Vv.op® ’Vo/^ 


(A- 37) 


Since the vapor partial pressure can be calculated from the gas lav 


Py,- ' Vv 


(A- 38) 


and the vapor peortial pressxure ratio, Y^, is 


JLo 

p_ 


(A-39) 


it follows that: 


'l - \ 

3 


ii. 

_Y^V J 

“ 30 

JPy,o 



(A-J»0) 


The Fourier ti^nsforms of Eqs. (A-26), (A-ll), (A-32), and (A-l»0) cwi 
be written in matrix form as 
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(A-l»2) 


^ is the system state vector and tj is the syston forcing function. 
The solution is then 

f j * j/‘^U 


(A-l*3) 
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After solving (A-i»3)« the entropy source trsnsfer function is 
found using Eq. (A-IO) 


. ^ . 



■^TL 




jUmI ^ 

p 


YP, 


(it-W) 


and the 
Syahols 

si 

\ 

°P 

®AB 

H 

i 

(MW) 

®L 

n 

P 

Pr 

Re 

r 

JM 


Bass source transfer Tju rtion is found using Eq,. (2.1.8). 
Used in This Appendix 
matrix 

specific heat of liqxiid fuel 
gas specific heat at constant pressure 
Bass diffusivity of A through B 
diameter, m 
enthalpy 
( -1 

molecular weight 

liquid fuel droplet mass, kg 

heat transfer Nusselt number, hd/»Cj^ 

mass trMisfer Husselt nijmber, kg/D^ 

manber of peu*ticles per unit volume 

pressure, N/m 

Prandtl number, ucp/tc 

gas constant 

Reynolds ntnnber, du/u 

radius, m 

transfer function for response of entropy source to a pres- 
sure perturbation 


Sc 


Schmidt number, w/D 


AB 


•<<f M C C 4 
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• Mitroi^ of e^* J/kg^K 

* 

t teB3>«rftture» K 

lyttea forcing function 
ipeedt m/sec 
Cartes i Ml coordinate 
system state vector 
¥ relative concentration « partial pressure ratio or density 

ratio 

molecular weight ratio of vapor to gas 

ratio of specific heat at constant pressure to specific heat 
at constant voIubk 
latent heat parameter, Wp 
gas thermal conductivity 
liquid droplet suiss fraction, 
tism, sec 
gas viscosity 

3 

density, kg/m 

collision diameter a Lennard-Jones parameter, A 
diffusion relaxation time, sec 
liquid droplet thermal relaxation time, sec 
^ mass source rate 

collision integral based on Lennard-Jones Potential 

Iv 

u angular frequency, radians/sec 

Subscripts: 

D diffusion 

H heat transfer 


A,B 


TL 


fuel 'vapor 

fuel r^r - liquid fuel 
ftwl YBp - air 
■een value 
perturbed value 

far froB soot x>artlcle or droplet 
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SPECTRUM FUNCTION CONFIDENCE UMITS 
purpose of this appendix is to discuss 8<aae of the statisti- 
cal and ec^putational considerations involved in the random (teta analy- 
sis reported herein. The primary source for this material is Bendat 
and Piersol (1971). 

The digital signal processor performing the random data analysis 
discussed herein digitised each saiiq>le record of data into 103lt words 
at a sampling rate* f , 2,oU8 times the highest frequency selected for 
analysis » f^. Accordingly* the processor mea>ory period or sample 
record length is 

T « (3m)2SLS.t / i \ , Ipgtt , ^ sec 

V Record / \Nunber of words/sec/ 2.0^8 f f * Record 

The resolution bandwidth Is 

, r 

g * i * — M. 
e T ‘>00 

Next, the finite Fourier transform is taken to obtain a single es- 
timate of the power spectral density function 

\(t) - I |x(r,T)i^ 


120 


where 
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X(f.T) - f dS - Xj,(f,T) ♦ iXj(f.T) 


and 


|x(f,T)(® » xJ(f,T) ♦ xJ(f.T) 

The sample distribution of G^(f«T) is detensined as follows. 
Assiadog X(t) is a Gaussian randc» variable, then Xj^(f,T) and 
Xj(f,T) are Gaussian uncorrelated randc» variables with zero neam and 
equal variances since the Fourier transfotm is a linear operation. 
Conseqxiently, the quantity 


|x(f,t)l^ * x|(f,T) + Xj(f,T) 


Is the suBi of the square of tvo independent Guassian variables. Thus, 
each frequency component of the estimate G(f,T) will have a sampling 
distribution 


0441 

G(f) 



vhere the chi-square variable with n = 2 degrees of freedom. 

The (l - a) confidence interval for the single record is 


x|(a/2) 


< G(f) < 


2 

X-, 



The (1 - a) confidence interva.1 expressed in terms of decibels is 
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10 logj^QO(f) ♦ 10 < 10 log^QG(f) < 10 logjQO(f) 


+ 10 


m 

- .«)] 


For a 9^ percent confidence interval 


x|(a/2) » X (0.005) = 10.597 


X^d - a/?) = X^(0.995)0.01 


Thus 


j^lO log^QC(f) + 10 log^Q ^3^ < 10 log^oG(f) < 10 log^QG(f) 


10 log. 


10 


O.OlJ 


or 


[lO log^QG(f) - 7.2 dB < 10 log^QG(f) < 10 log^^GCf) + 23 dfij 

'ilie error for a single estimate is unacceptable. The error is 
reduced by averaging over an ensemble of estimates. This is done by 
computing N disjoint, that is, independent, sample records and av^ 
eraging the N estimates at each frequency of spectral con 5 >onent. 

The total number of degrees of freedom becomes 


n = 2N 
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For a total tape record of 120 seconds » the Daxlnm ntaaber of indepen- 


dent saaple records is 


H « 120/T 


For the conputation reported herein: 


f = 1000 Hz 

B1 


f » 20U8 Hz 
s 


T = 1/2 sec 


= 2 Hz 


N » 2U0 


However, due to the processor design emd the use of redundant data the 
number of independent samples was less than 2l*0 but more than 128. For 
the following calculations N = 128 is used. Thus 


n = 256 


For n greater than 30, a normal approximation to the chi-square 
distribution is adequate. The expression 




'2n - 1 


is approximately normally distributed as the standard norm distribu- 


tion. Hence 


(z < X^(a)) » O = ♦ - 


For the 99“Percent confidence interval 


12U 


^2x^(o/2) - - I • 2 . 58 

y / 2x^(1 - o/2) - y2n - 1 » -2.58 


Thus 


X^(l - «/2) ■ [-2.58 + ^255]^/2 = 89.6289 

X^(o/ 2 ) - [2.58 ♦ ^ 255 p /2 » 172.0275 

Consequently « the 99“pe*‘cent confidence interval for this 128 point 
exaaple is 

[10 loei(, 0 (f) . 10 ‘ “ 


+ 10 log 


128 

89.(^289 


or 


1 10 log^pGCf) - 1.3 


dB < 


10 log^QG(f) 


10 log^pG(f) 


+ 1.5 dB 


BIBLIOGRAPHY 


Abdelhaaiidt A. N. , Harrje, D. T. , Plett, E. G. , and Sunanerf ield » M. » 
"Noise Characteristics of Combustion Augmented High-Speed Jets," 

AIAA Paper No. 73-189, 1973. 

Bendat, J. S. and Piersol, A. G. , lUindom Data ; Analysis and Mesisure- 
ment Procedures . Wiley-Interscience, 1971, p. 85. 

Bergman, P. G. , "The Wave Equation in a Medium With a Variable Index 
of Refraction,” J. Acoust. Soc. Am. . Vol. 17* No. 1*, April 19*»6, 
pp. 329-333. 

Bird, R. B. , Stewart, W. E. , and Lightfoot, E. N. , Transport Phenomena . 
John Wiley, New York, I960. 

Bonnell, Marshall, R. L. and Riecke, G. T. , "Combustion Instability in 
Turbojet and Turbo fan Augmentors," AIAA Paper No. 71-698, 1971. 

Byron, P. W. and Puller, R. W. , Mathematics of Classical and Quantum 
Physics . Addison-Wesley Publishing Co, , 1970, p. 395. 

Candel, S. M. , "Acoustic Conservation Principles and an Application to 
Plane and Modal Propagation in Nozzles and Diffusers," J. Sound 
Vib. . Vol. i*l. No. 2, 1975, PP- 207-232. 

Chow, J. C. F. , "Attenuation of Acoustic Waves in Dilute Emulsions and 
Suspensions," Journal of the Acoustical Society of America . Vol. 36 , 
December 196J<, pp. 2395-21*01. 

Cole III, J. E. and Dobbins, R. A., "Propagation of Sound Through At- 
mospheric Fog," Journal of the Atmospheric Sciences . Vol. 27, 1970, 
pp. l*2(S-l»3l*. 


125 


126 


Cole III» J. E. and Dobbins « R. A. , ^'MeasurcAents of the Attenuation 
of Sound in a Warn Air Fog,” Journal of the Atnospheric Sciences . 
Vol, 28, March 1971, pp. 202-209. 

Conte, S. D. and deBoor, C., Eleaant a ry Bvanerical Analysis , McGraw- 
Hill, Hew York, 1972, pp. 7^-83. 

(himnings. A., "Sound Generation and Transmission in Flow Ducts With 
Axial Temperature Gradients," J. Sound Vib. . Vol. 57, No. 2, 1978, 

pp. 262-279. 

Davidson, G. A., "Sound Propagation in Fogs,” Journal of the Atmo- 
spheric Sciences . Vol. 32, November 1975, PP. 2201-2205. 

Dobbins, R. A. and Tenkin, S. , "Propagation of Sound in a Gas Particle 
Mixture and Acoustic Combustion Instability,” AIAA Journal . Vol. 5» 

1967. pp. 2182-2186. 

Duff, G. F. and Naylor, D. , Differential Equations of Applied Mathe- 
matics . John Wiley & Son, Inc., 1966, pp. 

Qnmerling, J. J. , "Experimental Clean Combustor Program, Phase I, 

Noise Measurement Addendum," NASA CR-13i<853 (July 1975). 

Bmmerling, J. J. and Bekofske, K. L. , "Experimental Clean Combustor 
Program, Phase II, Noise Measuranent Addendum," NASA CR-1350l<5 
(January 1976). 

Epstein, P. S. and Carhart, R. R. , "The Absorption of Soxind in Suspen- 
sions and ftnulsions. I. Water Fog in Air," Journal of the Acous- 
tical Society of America . Vol. 25, May 1953, pp. 553-565. 

Eversman, W. , "Acoustic Energy in Duct, Further Observations," J . Sound 
Vib., Vol. 62, No. I4, 1979, pp. 571-532. 


127 


Faetht G. M. , "Current Status of Droplet and Liquid Coabuation," Paper 
presented at the Coabustion Institute. Spring Technical Meeting. 
Cleveland. Ohio. March 28-^. 1977« 

Friedaan. B. , Principles and Techniqties of Applied MathM»tlc8 « Wiley. 
Hew York. 1956. pp. l6U-l67. 

Goldstein. M. E. , Aeroacoustlcs . IfcGraw-Hill International Book Co., 
1976. pp. 30. 263. 

Gouldin. P. C. . Coabust. Sci. Technol ., Vol. 7. 1973. pp. 33-^*5* 

Graves. C. C. and Bahr. D. W. . "Atomization and Evaporation of Liquid 
Fuels." in Basic Considerations In the Combustion of Hydrocea*bon 
Fuels With Air , ed. Barnett, H. C. and Hibbard. R. R. , HACA Report 
1300. 1959. pp. 1-31. 

Harrje. D., ed. , "Liquid Propellant Rocket Combustion Instability," 

NASA SP-19**. 1972. 

Hurle, I. R., Price, R. B., Sugden, F. R. S., T. M. , and Thomas, A., 
"Sound Etaiosion from Open Turbulent Flames," Proc. Roy. Soc. A , 

Vol. 303, 1968, pp. U09-l»27. 

Igeurashi, J. and Toyama, M. , "Fundamentals of Acoustical Silencers. 

I - Theory and Experiment of Acoustic Low-Pass Filters," Rep. No. 339 
Aeronaut. Res. Inst., Univ. of Tokyo, Vol. 2 I 4 , No. 10, 1950. 

Ingard, U. and Singhal, V. K., "Effect of Flow on the Acoustic Reso- 
nances of an Open-Ended Duct," J. Acoust. Soc. Am. . Vol. 58 . No. U, 
October 1975, pp. 708-793. 

Jones, Robert E., "II. Emissions Reduction Technology Program," in Air- 
craft Engine Emissions . NASA CP-2021, 1977, PP. 22-23. 



128 


KaretaMr* A. and Seahotko, N. « "Core Molsa Source Dlagnoatlca on a Tur» 
bofan Engine Using Correlation and Coherence Techniques »** NASA 
TO X-73535. 1976. 

Karchsttr, A. N. , Resbotko, N. « and Nontegani* F. J.« **Nea8ureaent of 
Far Field &»bustion Eoise From a Turbofan Engine Using Coherence 
Functions," AIAA Paper No. 77-1277, 1977. 

Karchner, A. N., "Identification and Measurenent of Combustion Noise 
From a Turbo fan Engine Using Correlation and Coherence Techniques," 
Ph.D. Thesis, Case Western Reserve Univ. , Cleveland, Ohio, Jan- 
uary 1978. 

Kasin, S. B. and Enaoerling, J. J., "Low-Frequency Core Engine Noise," 
ASME Paper 7*»-WA/AER0-2, 1971*. 

Knott, P. R. , "Noise Generated by Turbulent Non-Premixed Flames," AIAA 
Paper No. 71-732, 1971. 

Kraft, R. E. and Wells, W. R., "Adjointness Properties for Differen- 
tial Systems With Eigenvalue-Dependent Boundary Conditions, With 
Application to Fluw-Duct Acoustics," J. Acoust. Soc. Am. , Vol. 6l, 
No. k, April 1977, pp. 913-922- 

LaBq>ton, M. , "Transsiission Matrices in Electroacoustics," Acoustlca . 
Vol. 39. 1979, pp. 239-251. 

Linden, L. H. and Heywood, J. B., "SnK>ke Emission From Jet Engines," 
Combust. Sci. Technol. . Vol. 2, 1971, PP. 1»01-1»11. 

Lumsdaine, E. and Ragab, S., "Effect of Flow on Qmsi-One-Dimensional 
Acoustic Weive Propagation in a Variable Area Duct of Finite Size," 

J. Sound. Vib. . Vol. 53, No. 1, 1977, pp. U7-61. 

Lumsdaine, E. , "Calculation of Pressure Reflection Ratio," J. Sound 
Vib. , Vol. 52, No. 1, 1977, pp. U»5-ll*7. 


129 


Haling, 0. C. , **SiBqplified Analysis of Rijke PhenosMnon," Journal of 
the Acoustical P c L ^ ty of Ansriea , Vol. 35, 1963, p. 1058. 

Mani, R., "Discrete /Tequency Noise Generation From an Axial Flow Fan 
Blade Row," ASMS Trans.. Journal of Basic Engineering . Vol. 92, 

1970, pp. 37-»»3. 

Marble, F. E.. "S(»e Gasdynaaic Problems in the Flow of Condensing 
Vapor," Astronautics Acta . Vol. lU, 1969, PP* 585-6lU. 

Marble, F. E. and Wooten, D. C. , "Sound Attenuation in a Condensing 
Vapor," The Physics of Fluids . Vol. 13, No. 11, Novwaiber 1970, 
pp. 2657-266U. 

Marble, F. E. and Candel, S. M. , "Acoustic Attenuation in Fans and 
Ducts by Vaporization of Liquid Droplets," AIAA Journal , Vol. 13, 
Ho. 5, May 1975, pp- 63l»-639- 

Mathews, D. C. and Rekos, Jr., N. F. , "Direct Combustion Generated 
Noise in Turbopropulsion Systems Prediction and Measurement," AIAA 
Paper Ho. 76-579, 1976. 

Mechel, von F. , Schilz, W. , and Dietz, J., "Akustische Impedanz Einer 
Luftdurchstromten Offhug," Acustica, Vol. 15, 1965, PP- 199-206. 

Miwa, T. and Igarashi, J. , "Fundamentals of Acoustical Silencers. 

II - Determination of Four Terminal Constants of Aco\istical Ele- 
ments," Rep. No. 3^^. Aeronaut. Res. Inst., Univ. of Tokyo, 

Vol. 25, No. 1», 1959. 

Hohring, W. , "Energy Flux in Duct Flow," J. Souna Vib. . Vol. li'. 

No. 1. 1971, pp. 101-109. 

Morfey, C. L., "Sound Transnission and Generation in Ductb vith Flow," 
J. Sound Vib. . Vol. lU, Ho. 1, 1971. pp. 37-55. 


130 


Norfagr* C. L. • "Acoustic Energy in Ron-Ttaifbra Hovs," J. Bound Vib. « 
Vol. Ikt Ho. 2, 1971* PP> 139-170. 

Mornt.., M. 1'. luiii F<-:jhbutjk, II. » Itethodo of Theoretical Physics . KcOmw- 

Hill Book Conpsny« 1953* p. 793. 

Mungur* P. snd Pluablee* Jr.» H. E. » "Influence of the Jet Exhaust Flow 
Field on the Acoustic Radiation Impedance of a Jet Pipe Opening," 
AIAA Paper Eo. 79-0676, March 1979* 

Muthukrishnan, M. , Strahle, V. C., and Handley, J. C. , "Effect of Coni- 
cal Flameholders on Combustion-Generated Noise," AIAA Journal, 

Vol. ll». No. 8, 1976, pp. 995-996. 

Norgren, C. T. , "Determination of Primary-Zone Smoke Concentrations 
From Spectral Radiance Measvirements in Gas Turbine Combustors," 

NASA W D-61*10, 1971. 

Parrott, Tony L., "An la^roved Method for Design of Expansion Chaniber 
Mufflers With Application to an Operational Helicopter," NASA 
TN D-7309, 1973. 

Plett, E. G., Leshner, M. D. , and Summerfield, M. , "Combustion Inten- 
sity and Distribution Relation to Noise Generation," AIAA Paper 
No. 75-5214, 1975. 

Polis, M. P. and Goodsen, R. E. , "Parameter Identification in Distribu- 
ted Systems: A Synthesizing Overview," Proceedings of the IEEE . 

Vol. 61*, No. 1, January 1976, pp. l»5-6l. 

Reshotko, M. and Karchmer. A., "Combustor Fluctuating Pressure Measure- 
Mcts In-Engine and in a Component Teat Facility - A Preliminary 
Comparison," NASA TM-738U5, 1977. 


ReshotkOt N. » Ktrcli— r» A. • Prako* P. F. * and NeArdle* J. 0«t Core 
Roita NiMur«Matt on m Yf^lOP Tuz^lte AIAA P^r Ro. 

77-21 I 1977. 

8hiaMlMak«r«» B. R.* ”0 m Turbina Engine Gore Roisa Source Isolnticm 
by latamal-to-Fnrfiald Correlations," AIAA Paper Ro. 77-1276, 

1977. 

Siaons, 0. A. and Lewis, P. F. , "Mms Transport and Heterogeneous 
Reactions in Porous Nediuai," Paper presented at the Coabustion 
Institute, Spring Technical Meeting, Cleveland, Ohio, March 2^30, 
1977. 

Saith, T. J. B. and Kilhaa, J. K. , "Roise Generation liy Open Turbulent 
Flaaes," The Journal of the Acoustical Society of Aaerica , Vol. 35* 
Bo. 5. m 1963, PP. 71V72H. 

Sofrin, T. G. and Ross, D. A. , "Experiaental Clean C(»nbU8tor Program, 
Phase I, Roise Addendua," RASA CR>13U820 (October 1975). 

Sofrin, T. G. and Riloff, Jr., R., "Experiaental Clean Combustor 
Program - Roise Study," RASA CR-135106 (Septeaber 1976). 

Stevart, W. L. and Weber, R. J., "A Review of NASA's Propulsion Pro- 
grams for Civil Aviation," AIAA Paper Ko. 78-l»3, 1978. 

Stone, J. R., "On the Effects of Flight on Jet Engine Exhaust Boise," 
RASA TN X- 71619, 1975. 

Strahle, W. C. , Muthiikrishnan, M. , and Beale, D. H., "Coherence Be- 
twe«i Internal and External Noise Oeneratod by Gas Turbine Cc«bus- 
tors," AIAA Paper No. 77-20, 1977. 

Strahle, W. C., Mathukrishnon, M. , and Neale, D. K.. "Experiaental and 
Analytical Separation of Hydrodynamic, Entropy and Direct Combustion 
Roise in a Gas Turbine Combustor," AIAA Paper Ho. 77-1275, 1977. 



132 


Svin^Mskti N. A* * ”n» tewd Fltld (teiMnt«d tor ft Sowe* XUstritoutloii 
la ft Lone Duet Ounrliig Stewtd Flov,” J« Somd Ylb«. Vol. )»0. 

■o. 1» 1975* ». 51-76. 

Iftkfthftfthit 7. • &3>binft* N. J. * ftad Atmlaader* D. W. . Control ftnd Bramlfc 
Syi^w. AddisOffWesl^ Publifthing Co.« 1972» p. 266. 

Toftirliii S. and Oobbins* R. A., "Atteaufttiem and Dlaperaicm bgr 

Particulate-Relaxation ProceaBes.** Journal of the Acouatlcal Society 
of Aaerica. Vol. bO, Ho. 2, 1966. pp. 317-32U. 

Tetter. B. J.. "Acouatic Energy Flov in Lined DiKts Onatalning Utalfom 
Flou," J. Sound Yib. . Vol. 28. Ho. 2. 1973, pp. 205-21?. 

Ton Olatan. U. B. . "Correlation of Ooabustor Acouatie Power Lewela In- 
ferred Froaa Internal nuetuating Preasure Neaaiueaen^ " HASA 
HI-78986. 1978. 

Walty. J. R.« Vicln. C. E., and Wilaon, R. E. . »»n^ ala of Moaen- 

tum- Heat and Mass Trnnafer . John Wiley. Hew York. 1S>69. 

Vri^t. P.. "The Oxidation of Soot by 0 Atone." International Swapoalua 
on Coaibuation. 15th (Conbuation Inatitute. Pittsburg. Pa.. 1975). 
pp. lUk9-lh59* 


